{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "5c4a8f18",
   "metadata": {},
   "source": [
    "### 相位响应\n",
    "#### MATLAB® 函数可用于提取滤波器的相位响应。在给定频率响应的情况下，函数 abs 返回幅值，angle 返回以弧度为单位的相位角。要使用 fvtool 查看 Butterworth 滤波器的幅值和相位，请使用："
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 270,
   "id": "da6728f8",
   "metadata": {},
   "outputs": [],
   "source": [
    "from scipy import signal\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np             "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 271,
   "id": "f081f93c",
   "metadata": {},
   "outputs": [],
   "source": [
    "b, a = signal.iirfilter(9, Wn = 400,btype='lowpass', analog=True, ftype='butter', output='ba')\n",
    "w, h = signal.freqs(b, a,worN=np.linspace(0,1000,10000))    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 272,
   "id": "e770ca82",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, 'Phase(radians)')"
      ]
     },
     "execution_count": 272,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbEAAAEWCAYAAADoyannAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABMQUlEQVR4nO3dd3gUVffA8e8JoUuRJkWKNEWKdBVEkKJIkWKhiIAVFMWG7ZXXGP35ymtvrwIqgiJNpEgTRESQakSqdBFBQYpIb0nO74870TWmkWyY3c35PM882Z3ZO3Nms9mTe+fOvaKqGGOMMeEoyu8AjDHGmMyyJGaMMSZsWRIzxhgTtiyJGWOMCVuWxIwxxoQtS2LGGGPCliWxCCAi60SkRRD395OItA7W/iKFiFwsInFpbB8pIv+XwX3lFZENIlIqeBGmebxKIqIiEp2JsvNF5I7siMuYrLIklkXeF/4pESmRbP1K70ujUnbHoKo1VXW+d9ynRWR0dh8zh3oWeCkjLxSRFiKSKCJHvOUXEYlN2q6qJ4ERwGPZFOsZ8T7Hx71YfxORD0TkHL/jSuL9LR0NeC9fEZFcfsdl/GdJLDi2AT2SnohIbSC/f+GEhsz81x+qRKQMcBUw5QyK/aqq56jqOcAVwO0i0jlg+xigj4jkDVqgWdPRi7U+0AgY7HM8yV3ixdcc6Abc5nM8JgRYEguOj4DeAc/7AB8GvkBE2ovI9yJySER2iMjTybb3FpHtIrJfRP4d2KTn1a4miMiHInLYaz5sGFD2JxFpLSJtgX8B3bz/WFcFbg94/d9qayJyS8Cxn0wWV5SIPC4iW73tE0SkWEpvglf72Ckij4nIbuCDtMqLSD4RGe2t/0NEvhWR87xt80XkeRFZLiIHRWRq4HFF5DrvffjDe22NZO/HIBFZ7ZUdLyL5vG0lRGS6V+53EVkoIlHetrIi8qmI7BWRbSIyMOD02gArVPVEwHHqicgK73cyHsiX0vsCoKrbgMXAxQHrdgIHgMtSeT8bi8gSL9ZdIvKWiOQJ2K4i0l9ENovIARH5n4iIty2XiLwkIvtE5EegfWqxpRDrL8AsoFbA6ooissg71zkS0PIgIp+IyG7vvV4gIjUDtrUTkR+8cr+IyKCAbR3EtVj8ISKLRaROBuPbAiwC6mZkX97n8Rcvho0i0spb/7SITPQ+H4e93+UlAeVqeJ+tP7zP2nUB20Z67/cMr+wyEanibRMReVVE9njvyWoRqeVty+v9Xn4WV+MdKiI5/h/eLFFVW7KwAD8BrYGNQA0gF7ADqAgoUMl7XQugNu4fhzrAb0Bnb9vFwBHcf+t5cE1Wp4HW3vangRNAO2//zwNLk8cQ8NrRKcUY8PzP1wQc+0ogL/AKEB+wvweApcD53vZhwNhU3osWXtn/eq/Nn1Z5oB8wDSjgnVcDoLC3bT7wC+6LtCDwaUDM1YGjuMSSG3gU2ALkCTjf5UBZoBiwHujvbXseGOqVyw00A8T7vXwHPOX9DioDPwLXeOVeBP4XcK55gO3Ag95+bvB+Z/8X8F7sDHh9Ne98WiZ7zz4DBqbyfjbAJbhooJJ3Hg8EbFdgOlAUqADsBdp62/oDG4Dy3nvwlff66LQ+x97j8sA64NmA38VW733P7z0fElD2NqCQ9/t9DVgZsG0X0Mx7fC5Q33tcH9gDXOr97vt4MeRNJT4FqnqPL/L2+2B6+wIuxP09lvVeWwmoEvB3cNr73eUGBuFaVZI+G1tw/xTmAVoCh4ELvbIjgd+Bxt7v52NgnLftGtxnqSjus1UDKONte837nRfz3rNpwPN+f4+F8+J7AOG+8FcSG4z7gmwLfOF9sP9MYimUew141Xv8FAGJAfelfoq/J6a5AdsvBo4njyHgtWeSxJ5K+uPznhdMduz1QKuA7WW8P/x/fBnivrhPAfkC1qVaHvfltxiok8K+5vP3L8qLvX3nAv4NTAjYFoVLEC0CzrdXwPYXgKHe42eAqXhfiAGvuRT4Odm6J4APvMfvJovnSuBXQALWLebvSSwR+AM45H0WJuEl2oAyHwNPZfCz9gAwOeC5AlcEPJ8APO49noeXuL3nV5N+EjvixbsdeBvIH/C7GBzw2nuAz1PZT1HvOEW85z/j/lkpnOx17+AlyYB1G4HmqexXvffxqPd4LF7CS2tfQFVcgmsN5E72mqf5+z+DUXhJ11t2A1EB28cCT3uPRwLvBWxrB2zwHrcENuH+AQksL178VQLWXQ5sy8jv35aUF2tODJ6PgJ5AX5I1JQKIyKUi8pXXVHUQ959yUpNMWdx/iwCo6jFgf7Jd7A54fAzIJ8G55pT82EeTHbsiMNlrUvkDl5QSgPNS2d9eDWhyS6f8R8BsYJyI/CoiL4hI7oCyOwIeb8f9d1zCi3l7QMyJ3mvLBbw++fuV1EnhRdx/2HNE5EcReTwgzrJJcXqx/ivgPA/g/nNOUhb4Rb1vooAYA/2qqkVVtTDuy/04MCrZawrhEsc/iEh1r+lzt4gcAv7DX5+Z9M7zb7/XFGJLSWcv3oqqeo+qHk/vOF6z5RBxzcWHcMmQgDivx33BbxeRr0Xkcm99ReDhZO93eS/u1NT3jtsN909HwfT2pa7p8QFcwtojIuNEJPAYgZ/9RGCnF0NZYIe3Lsl2MvAZU9V5wFvA/4DfRGS4iBQGSuL+Qf0uIM7PvfUmkyyJBYmqbsc1RbTD/ced3BhcM0J5VS2Ca9ISb9suXHMbAF4befHMhpLCuqO4P54kpQMe78L9wScdu0CyY+8ArvW+3JKWfOqum2Tk+KmWV9XTqhqrqhcDTYAO/P3aYvmAxxVwNbh9uBpQxYCYxXttajH9FZzqYVV9WFUrAx2Bh7xrJDtw/xEHxllIVdt5RVfjmtOS7ALKJV2DCogxteMexH0GOibbVANYlUqxd3BNgtW8RPgv/vrMpOdvv9e0YsuinkAnXE2nCK65Drw4VfVbVe0ElMJ1ipngbd8BPJfs/S6gqmPTOpg6E4AluFaEdPelqmNU9Qr+auL/b8AuAz/7Ubi/w1+9pby3LkkFMvAZ8475hqo2AGriPjeP4D67x4GaAXEWUddZxWSSJbHguh13zeNoCtsKAb+r6gkRaYz7408yEegoIk28C/exZPzLKrnfgErJ/vhWAt1FJLe4DiE3JDt2BxG5wjv2M/z9czEUeE5EKgKISEkR6XQG8aRaXkSuEpHa4rpKH8IlqYSAsr3E3ZtVwItroqom4L4I24tIK6/m9jBwEteclyavA0BVL/kc8o6XgLuGdsjrBJDfq2HUEpFGXtEvgPridRDBfYnGAwNFJFpEuuKuj6R23HOA7rhrTUnryuGujSxNpVghL8YjInIRcHd65xdgghfb+SJyLvB4egUyqRDuvd+P+0fpP0kbRCSPiNwsIkVU9TR/vd/gmmf7ey0UIiIFxXV+KpT8AKkYAtwlIqXT2peIXCgiLcX1AD2BSyKBn7EGItLVa9V4wDuXpcAy3D9/j3p/Ny1w/4CMSy8wEWnkxZLb28cJIMGr1b0LvCre/YEiUk5ErsngOZsUWBILIlXdqqqp3Qx7D/CMiBzG/Qc5IaDcOuA+3B/ILtwF5D24P6gz9Yn3c7+IrPAe/xuogmsSi8XVCAKPPcBbt8t7zc6A/b2Oq0HO8WJfimvKyai0ypfGJdFDuGbGr4HAe9w+wl172I3r+TfQi3kj0At4E/ffbUdc9/BTGYinGjAXd/1nCfC2qs73kmNHXI+3bd5+38PVLlDV33DXmTp5z08BXXHNxwdwTVzJa+BlxbtPDNcUVQy4OWB7T2CUunvGUjLIe81h3Jff+AycX5J3cU21q4AVKcQWLB/izu0X4Af+mZBvAX7ymhr7435veH8nd+Ka3Q7gmnj7ZvSgqroG93l5JJ195cUlvH24z1EpXI02yVTc7+6AF2tXr4XgFHAdcK1X9m2gt6puyEB4hXHv/wHce7Ofv+4vfMyLb6n3nszFdT4xmSR/b9I3ocD7r/0PXDPSNp/D8YWIzMd1PnnP71iSiMjFuGtajTWLfzhezWAVcKWq7glGfObMiLvNpaqq9vI7FpN5VhMLESLSUUQKiEhB3H9ta/jrIrkJAar6g6o2ymoC8/Z1UlUvsgRmTNZYEgsdnfjrgnI1oHswviyNMSaSWXOiMcaYsGU1MWOMMWErIgZojYqK0vz5bfgxY4w5E8eOHVNVDevKTEQksfz583P0aEq3ZhljjEmNiBxP/1WhLawzsDHGmJzNkpgxxpiwZUnMGGNM2LIkZowxJmxZEjPGGBO2QjaJiUhbcVOJb5G/5nwyxhhzFoX6d3FIJjFvao7/4UaQvhjo4Q2+aowx5iwJh+/iUL1PrDGwRVV/BBCRcbixBX8I5kG2bD/OzSMep5n+mwLeRLSBUxxm5HGolQm1eFIrExUF+fNDgQL/XIoVgxIlIFcujDH+OivfxVkRqkmsHH+fWn0nyeawEpG7gLsA8uTJk6mDzNsQx/LEYSw/MhXGT4Jd9TMZrgm2qCiXyEqVgkqV4KKL3FK7NtSrB7lz+x2hMREhWkQC50AcrqrDA56n+13st1BNYinNavy3kYq9N3o4QMGCBTM1ivFd1zSjXq2FdJ3QlX0DmvJex/e4uc7N/3hd4BjJycdLTm1buJTxK56EBDhxAo4d+/ty9Cjs3w+//fbXsm0bfPEFnPSmjixQAC67DK65Bm68ES64AGNM5sSrasM0tqf7Xey3UE1iO4HyAc/Px01REnSNyjXiu7u+48ZPbqTX5F7E/RrHi1e/SHTUX29NWs1j5uxISIDt2+G772DhQliwAB57zC2NG8N998FNN0EmK+XGmJSdte/izArJqVhEJBrYBLTCTXv+LdBTVdel9PqCBQtqVsdOPJ1wmkFzBvHG8jdoUakFE26YQMmCJbO0T5O9tm2DiRPhgw9g/XooXRr+/W+4805rbjQmI0TkmKoWTGP7GX0X+yEkkxiAiLQDXgNyASNU9bnUXhuMJJbkw1Uf0m96P0oWKMnkbpNpULZBUPZrso+qa2587jlXQ6teHd55B1q29DsyY0JbeknMe02Gv4v9ELJJ7EwEM4kBfPfrd3QZ34U9R/cwvONwel/SO2j7NtlHFaZPh4cegi1b4J574IUXoGCaf6LG5FwZSWKhLiTvE/Nbg7IN+O6u77i8/OX0mdKH+2fdz+mE036HZdIhAh07wqpV8OCDrjZ22WWwdavfkRljsovVxNIQnxjPo188yqtLX+XKilfyyY2fUKpgqaAfx2SPL76A7t0hMREmTIA2bfyOyJjQYjWxCBcdFc0r17zCR10+Yvkvy2kwvAHf/vKt32GZDGrTBuLioEIFaN8exo/3OyJjTLBZEsuAXnV6sei2ReSSXDT7oBkjV470OySTQRdc4Dp7XHYZ9OgBw4b5HZExJpgsiWVQ/TL1ibsrjqYVmnLr1Fu5b+Z9dp0sTBQpArNnQ7t20L8/jBjhd0TGmGCxJHYGShQowexes3nosod469u3aPVhK3478pvfYZkMyJ8fPv3UjfJx553uGpkxJvxZx45MGrNmDHd8dgfF8hdjUrdJNC7X+Kwe32TOsWPQti0sWQIzZsDVV/sdkTH+sY4dOVjP2j1ZfPticufKTbMPmjHie2ujCgcFCrh7yWrWdOMurl3rd0TGmKywJJYFdUvXJe7OOK6seCW3f3Y798y4h1MJp/wOy6SjcGGYNs3dBN2hgxtk2BgTniyJZVHxAsWZdfMsHmnyCO/EvUPLUS3ZfWS332GZdJQv7xLZnj3QqZMbUd8YE34siQVBdFQ0L7R5gbHXj2XFrhU0GN6ApTuX+h2WSUeDBjB6NCxbBgMH+h2NMSYzLIkFUfda3Vly+xLy5spL85HNeW/Fe36HZNLRtSv861/w7rvW9d6YcGS9E7PB78d/p8enPZizdQ79GvTj9bavkzc6r99hmVQkJLgeiwsXwuLFUN8m+DY5RCT0TrQklk0SEhN4ct6T/HfRf7n8/MuZeNNEyhYq63dYJhV797rmxagoN/Fm8eJ+R2RM9ouEJGbNidkkV1QuhrQewvgbxrPqt1U0GN6AxTsW+x2WSUXJkm6CzV274NZb3bQuxpjQZ0ksm91U8yaW3r6UgrkL0mJkC4bF2eB9oapxY/jvf12vxf/9z+9ojDEZYc2JZ8mB4wfo8WkPZm+dzV317+LNdm+SJ1cev8Myyai6OcnmznW9Fi+5xO+IjMk+kdCcaEnsLEpITGDwvMEMWTSEZhWaMfGmiTY/WQjau9clr8KF3fUxmxnaRKpISGLWnHgW5YrKxfOtn2dM1zF8++u3NHq3ESt3r/Q7LJNMyZLw8cewaRPcf7/f0Rhj0uJLEhORG0VknYgkikjDZNueEJEtIrJRRK7xI77s1qN2DxbeupCExASajmjKJ+s+8Tskk8xVV7n7x95/3ybTNCaU+dKcKCI1gERgGDBIVeO89RcDY4HGQFlgLlBdVRPS2l+4NCcmt+vwLrpO6MrSnUt56sqniGkRQ5RY5ThUnD4NzZvDunWwejVUrOh3RMYElzUnZpKqrlfVjSls6gSMU9WTqroN2IJLaBGpTKEyzO8zn751+/LMgme4YcINHDl1xO+wjCd3btesqAp9+riboo0xoSXU/u0vB+wIeL7TW/cPInKXiMSJSFx8fPxZCS475I3Oy4jrRvDqNa8ydeNUmrzfhG0HtvkdlvFccAG8+SZ8/TW8/LLf0Rhjksu2JCYic0VkbQpLp7SKpbAuxfZOVR2uqg1VtWF0dHRwgvaJiPDAZQ8w6+ZZ7Di0g0bvNuKrbV/5HZbx9O4NN9wAgwfD99/7HY0xJlC2JTFVba2qtVJYpqZRbCdQPuD5+cCv2RVjqLm6ytUsv2M5JQuWpM1HbXj727eJhFsgwp0IDB3qei3efDMcP+53RMaYJKHWnPgZ0F1E8orIBUA1YLnPMZ1V1YpXY+ntS2lbtS0DZg6g//T+NtFmCCheHEaOhPXr4bHH/I7GGP+JyIsiskFEVovIZBEp6kccfnWx7yIiO4HLgRkiMhtAVdcBE4AfgM+BAen1TIxERfIVYWr3qTzW9DGGrxhO6w9bs/foXr/DyvHatHH3jb35Jsye7Xc0xvjuC6CWqtYBNgFP+BGEjdgR4sasGcPtn93OeQXPY2r3qVxS2sZB8tPx49CoEezfD2vWQIkSfkdkTOYFq4u9iHQBblDVm4MQ1hkJteZEk0zP2j1ZeOtC4hPjaTKiCVM2TPE7pBwtf37X7f7336FfPxvt3oS96KRe3t5yVyb3cxswK5iBZZTVxMLErsO76DSuE3G/xjGk9RAeafIIIil15jRnw4svwqOPutmgb73V72iMyZz0amIiMhconcKmJ5M66YnIk0BDoKv6kFAsiYWR46eP03dqXyasm8CtdW9laIehNhK+TxISoHVriIuDlSuhShW/IzLmzGW1OVFE+gD9gVaqeix4kWWcNSeGkfy58zP2+rE8deVTfLDyA9p81IZ9x/b5HVaOlCsXjBrlft5yC4Tx/fbGZIqItAUeA67zK4GBJbGwEyVRxF4Vy+guo1m2cxmXvXcZG/Zt8DusHKlCBXjnHViyBJ5/3u9ojDnr3gIKAV+IyEoRGepHENacGMYW71hM53GdOZVwiok3TaR15dZ+h5Qj9ewJEybAokVw6aV+R2NMxkXCAMCWxMLcT3/8RIcxHdiwbwNvtXuL/g37+x1SjvPHH24Szdy53bBUhQr5HZExGRMJScyaE8NcpaKVWHz7Yq6ucjV3z7ibBz5/gITEHHd/uK+KFoXRo2HbNhg40O9ojMlZLIlFgMJ5C/NZj8+4/9L7eX3Z61w37joOnTzkd1g5SrNmbhLNkSNd06Ix5uyw5sQIMzRuKPfOvJcaJWswrcc0KhWt5HdIOcbp0y6ZbdwIq1a5jh/GhDJrTjQhp3/D/nze63N2HtrJZe9dRtyvcX6HlGMkTaIZH++mb7FJNI3JfpbEIlDryq1ZfNti8ufOT/ORzZm2cZrfIeUYVar8NYnmCy/4HY0xkc+aEyPY7iO76TCmA9/v/p432r7BgMYD/A4pR1CFbt1g8mRYvNgNGGxMKIqE5kRLYhHu6Kmj9Pi0B9M2TWPQ5YP4b5v/EiVWAc9uBw5AnTpuwOAVK+Ccc/yOyJh/ioQkZt9mEa5gnoJM7jaZAY0G8NKSl+g2sRvHT9vUxNnt3HNdt/stW+CBB/yOxpjIZUksB8gVlYs3r32Tl9q8xMQfJtL6o9Y25uJZ0Lw5PP44vP8+fPqp39EYE5msOTGHmfjDRHpN6kX5IuWZdfMsqhar6ndIEe3UKWjaFLZuhdWr4fzz/Y7ImL9Yc6IJOzdcfAPz+szjwPEDXP7+5SzZscTvkCJanjwwZgycPGnd7o3JDr4kMRF5UUQ2iMhqEZksIkUDtj0hIltEZKOIXONHfJGuSfkmLLl9CUXzFaXlhy2ZvH6y3yFFtGrVXLf7r76y0e6NCTa/amJfALVUtQ6wCXgCQEQuBroDNYG2wNsiksunGCNateLVWHL7EuqWrssNn9zA0DhfZlHIMW691Y12HxMDCxb4HY0xkcOXJKaqc1Q1aRrBpUDSlYJOwDhVPamq24AtQGM/YswJShQowZe9v+Taqtdy94y7iZ0fSyRcIw1FIjB0KFSu7JLZPutXY0xQhMI1sduAWd7jcsCOgG07vXX/ICJ3iUiciMTF27S6mVYgdwEmd5tMn0v68PTXT3PPjHtsFPxsUqiQGxx4717o2xcSE/2OyJjwl21JTETmisjaFJZOAa95EogHPk5alcKuUqwaqOpwVW2oqg2jo6ODfwI5SO5cufmg0wc83vRxhn43lJsm3sSJ+BN+hxWR6tWDl1+GGTPg1Vf9jsaY8Jdt3/6qmuY0wyLSB+gAtNK/2rB2AuUDXnY+8Gv2RGgCiQjPt36e8845jwdnP0jb0W2Z0n0KRfMV9Tu0iDNgAMyb5+4ha9YMGluDuTGZ5st9YiLSFngFaK6qewPW1wTG4K6DlQW+BKqpaprtW3afWHCNWTOGvlP6UqNkDWbdPIuyhcr6HVLEOXDA1cpE3GzQRYv6HZHJiew+scx7CygEfCEiK0VkKICqrgMmAD8AnwMD0ktgJvh61u7J9J7T2fr7VpqOaMqm/Zv8DininHsujB0LO3bAnXe6QYONMWfORuwwqYr7NY52H7dDUeb0mkO9MvX8DinivPACPPYYvP023H2339GYnCYSamKWxEyaNu3fRJuP2nDwxEFm9JxB0wpN/Q4poiQmQrt2MH8+LFnimhiNOVssiYUIS2LZ6+eDP9P6w9b8cvgXpnSbQpsqbfwOKaLs2QP160PevPDdd3Z9zJw9kZDEQuE+MRPiKhSpwMJbF1Ll3Cp0GNuBqRum+h1SRClVyt0/9vPP0KeP3T9mzJmwJGYy5LxzzmN+3/nUK12P6ydcz8erP06/kMmwJk3gxRfhs8/cT2PChYgMEhEVkRJB2V+sFJTYjA83aM2J5owcPnmYTuM6Mf+n+bzT/h36Neznd0gRQxW6dXNzj335JbRo4XdEJtJltTlRRMoD7wEXAQ1U9YwHVJNYicKNmXsz0Ag4CeQF9gIzgeEao5tTK59uTUyEUiJ0EWGACLeJ0FjEanA5VaG8hZjRcwbtqrWj/4z+vLjIqg3BIgLvvedGve/eHXbt8jsiY9L1KvAoqYyslEFfAVVwA8GX1hgtrzFaCmiGG1t3iMRKr9QKp1oTE+Eq4HGgGPA9sAfIB1T3DjgReFmVQ1kIPiisJnb2nU44zS2Tb2H8uvE83fxpYlrE+B1SxFi3zo3i0aCBq5Hlzu13RCZSicgpYE3AquGqOjyDZa/Djbh0v4j8BDTMZE0st8bo6cy+Jq0k9iLwpio/p7AtGjdkVC5VfJ943ZKYPxISE7hj2h2MXDmSp658iqdbPI1ISsNfmjP18cfQqxcMGmTXyEz2Sa85UUTmAqVT2PQk8C/galU9mJUk9uexYqUKsFNj9KTESgugDvChxugfaZaza2ImKxI1kTs/u5MRK0cwuNlgnrnqGUtkQXLPPfDOO+4aWdeufkdjIlFmr4mJSG3csIDHvFVJ49w2VtXdmYolVlYCDYFKwGzgM+BCjdF2aZVL89qWCM1FqOM9vkmEt0R4UIS8mQnSRJ4oieLd697ljnp38H8L/4/B8wbbnGRB8uqr0KiRm7Zlk438ZUKIqq5R1VKqWklVK+EGb6+f2QTmSdQYjQe6AK9pjD4IlEmvUKqj2IvwP1x1Lq8Im4BzcOMZNgFG4HqSGEOURDGs4zCiJIr/fPMfEjWR/7T6j9XIsihvXvjkE3dtrEsXWLrUzUlmTIQ6LbHSA+gDdPTWpXtFOK2pWK5S5WIR8gG/AKVUSRBhGLA6y+GaiBIlUbzT4R2iJIohi4aQqIkMaT3EElkWVawI48fD1Ve7G6EnToQo6xtsQoxXG8uqW4H+wHMao9skVi4ARqdXKK2OHStUqZ/8cUrP/WbXxEKHqnLvzHt5O+5tBl0+iBfavGCJLAheeQUefhiefRYGD/Y7GhMpImHYqbRqYqVEeAg323LSY7znJbM9MhOWRIS32r2FiPDSkpfInzs/z1z1jN9hhb0HH4QVK+Cpp9wgwe3b+x2RMcElsdIUeBqoiMtNAqjGaOW0yqWVxN7FzfmV/DG4O7SNSZGI8Ma1b3Ai/gTPLniWgrkL8tgVj/kdVlgTgeHD4YcfoGdPWL4cLrzQ76iMCar3gQeB74AMzyOZahJTJTYIQZkcKkqiGNZhGMdOH+PxLx+nQO4C3HfpfX6HFdYKFIDJk6FhQ+jcGZYtg8KF/Y7KmKA5qDE660wLpdU78Y20Cqoy8EwPZnKWXFG5GNV5FMdOH2Pg5wMpmKcgt9W7ze+wwlrFiq7HYuvW0Ls3TJpkHT1MxPhKYuVFYBJu/EQANEZXpFUorY//d96SD6gPbPaWupxBVc/kbLlz5Wb8DeO5pso13PHZHYxdM9bvkMJeixbw8sswdSr83//5HY0xQXMp7mbn/wAve8tL6RVKd8QOEb4CrlbltPc8NzBHlauyGnGwWO/E0Hfs9DGu/fhaFv28iEndJnHdhdf5HVJYU3U3QX/4oUtm19nbaTIhEnonZiSJbQQuV+V37/m5wFJVMn1ZWUSeBToBibiBhfuq6q/etieA23G1vYGqOju9/VkSCw+HTx6m1YetWLNnDXNvmUvTCk39DimsHT8OV14JGzbAokVQp47fEZlwE2pJTGKlPVAT1wIIgMZomt2bM9KaPgT4XoSRIowEVuCqe1nxoqrWUdW6wHTgKQARuRg3r0xNoC3wtkjGJ0czoS1pGpfyhcvTYWwH1u1Z53dIYS1/flcLK1wYOnaEPXv8jsiYzJNYGQp0A+7Dda+/EdfdPk3pJjFVPsC1VU72lstVGZWVYFU1cPqWgvw1F00nYJyqnlTVbcAWoHFWjmVCS8mCJZndazb5ovPR9uO27Di4w++QwlrZsi6R7d3rhqY6eTL9MsaEqCYao72BAxqjscDlQPn0CqWaxESolPRYld2qTPWW3d52EeH8zEYrIs+JyA7cGIxPeavLAYHfaju9dSmVv0tE4kQkLj4+PrNhGB9ccO4FfH7z5xw6eYhrRl/D78d/9zuksNawIYwaBYsXw113uetlxoSh497PYxIrZYHTwAXpFUqrJvaiCJ+K0FuEmt4MzxVEaCnCs8AioEZqhUVkroisTWHpBKCqT6pqeeBj4N6kYinsKsU/SVUdrqoNVbVhdHRa92ybUHRJ6UuY0m0KWw9spePYjhw7fSz9QiZVN94IzzzjOnq88ILf0RiTKdMlVooCL+IuW/0EjEuvUJodO0S4GFdTaoobEv8YsB6YCUxU5URWoxaRisAMVa3ldepAVZ/3ts0GnlbVJWntwzp2hK9P1n1Ct4nd6HxRZybeNJEosZueMkvVjeYxfry7KbpTJ78jMqEu1Dp2JJFYyQvk0xg9mO5r/Zj7SUSqqepm7/F9QHNVvUFEagJjcNfByuImXaumqmnel2ZJLLy9uuRVHprzEI83fZznWz/vdzhh7fhxaN7cDU+1aBFcconfEZlQFgpJTGKlpcboPImVFKd+1RidlFb5dNvhRCgAPARUUOUuEaoBF6oyPVMRO0NE5EJcF/vtuOH3UdV1IjIB+AGIBwakl8BM+HvgsgfYuH8jQxYNoXrx6txa71a/QwpbST0WGzVyPRaXL4fSKU0ub0zoaA7M4685xAIpbgSPVGXkPrHxuJE7eqtSS4T8wBJV6mYq3GxgNbHwdzrhNO3GtOPrn77mi1u+oHml5n6HFNZWrIArrnD3jn31lUtuxiQXCjWxrMpIEotTpaEI36tSz1u3SpWQaaiwJBYZ/jjxB5e/fzl7ju5h6e1LqVa8mt8hhbXJk+H666FrV5gwwcZYNP8UCklMYuWhtLZrjL6S1vaMfKxPebUvBRChCgGDMxoTLEXzFWV6j+kIQoexHfjjxB9+hxTWunRxYyx++ik8+qjf0RiTqkLe0hC4G3dbVTncZaaL0yuckZpYG2Cwt7M5uJ6KfVWZn5Wog8lqYpFl4faFtPywJW2rtmVq96nWYzELVGHgQHjrLbcMGOB3RCaUhEJNLInEyhzgeo3Rw97zQsAnGqNt0yqXkRE7vgC6An2BsUDDUEpgJvI0q9iM1655jembpvPs18/6HU5YE4HXXnOdPAYOhOlZ6Y5lTPaqAJwKeH4K/hp0IzWp1sREqJ9WQVXSnOPlbLKaWORRVW6deiujVo1iWo9pdKjewe+QwtrRo67r/fr1sGABNGjgd0QmFIRYTexJ4Cbc8IYKdAEmaIymOVZvWknsK+9hPlxb5SrciBp1gGWqXBGc0LPOklhkOn76OFd8cAVbf9/Kt3d+ax09smj3brjsMje+4tKlboJNk7OFUhIDkFhpAH/mlgUao9+nVybV5kRVrvLmDNsO1FeloSoNgHq4gXmNyVb5c+dn0k2TiI6Kpsv4Lhw5dcTvkMJa6dIwY4a7Ibp9e/jjD78jMubvNEa/w122mgzsl1ipkF6ZjFwxv0iVNX8eRFkLoXOPmIlsFYtWZNwN41i/bz13z7gbP0aYiSQ1a8KkSbBxo+t+f+pU+mWMORskVq6TWNkMbAO+9n7OSq9cRpLYehHeE6GFCM1FeBc3fqIxZ0Xryq2JaR7D6NWjGbUqS7MAGaBlS3jvPZg3D267DRIT/Y7IGACeBS4DNmmMXgC0xg00n6aMJLFbgXXA/cADuCGhbFwgc1Y92exJrqp0FQNmDuCHvT/4HU7Y69MHnnsOPv7Y7iEzmSci94nIRhFZJyJZnT/htMbofiBKYiVKY/QrMtDq58sAwMFmHTtyhl2Hd3HJ0Es475zzWHbHMgrkLuB3SGEt8B6yF1+EQYP8jsicbVnp2CEiVwFPAu1V9aSIlFLVTM8vLrEyF+gMPA+UAPYAjTRGm6RVLt2amAjbRPgx+ZLZQI3JrDKFyvBRl49Yu2ctD3z+gN/hhL2ke8huvBEeeQQ++sjviEyYuRsYoqonAbKSwDydcNN9PQh8Dmwl5UGB/yYjI3YUD3iaD7gRKKb652zMvrOaWM7yxNwnGLJoCGOvH0v3Wt39DifsnTwJ114LCxfCtGnQNs3xEUwkEZFT8FfHPWC4qg7PYNmVwFSgLXACGKSq32YqjljJBczWGG19xmUz05wowjd2n5jxy+mE0zQf2Zwf9v7AmrvXUL5Ieb9DCnsHD0KLFrB5sxv1vlEjvyMyZ0N6zYkiMhdIaTKfJ4HncFOo3A80AsYDlTWT16gkVj4DbsnIRJh/K5eBmljgyB1ReIM02ij2xk9bf99K3WF1aVS2EXN7z7XxFYNg925o0gQOH3YTalav7ndEJrtl8ZrY57jmxPne863AZaq6N1P7i5UJuN6JXwB/fqFrjA5Mq1xG/vJfDlieB+rjhgYxxjdVilXhtWte46ufvuK1pa/5HU5EKF0aZs9218quuQZ27fI7IhPipgAtAUSkOpAH2JeF/c0A/g0swM1hmbSkKSM1scqqf+/IIcIFqmzLfKzBZTWxnElV6Ty+M59v+Zy4O+OofV5tv0OKCHFxrmmxShWYPx/OPdfviEx2yWJNLA8wAtcN/hTumti8IIaXsTgykMRWqP59MGARvvOGoAoJlsRyrj1H91D7ndqcV/A8vr3zW/JG5/U7pIjwxRfQoYMbKPiLL6BgyIyuZ4IpFMZOlFiZBgwHPtcYPZ1sW2XcDCo/aYyOSKl8qs2JIlwkwvVAERG6Bix9cb0Usx68yCARUREpEbDuCRHZ4t1Ad00wjmMiV6mCpRhx3QjW7FnD4HmD/Q4nYrRpA2PHwrJlbnLNkzYNrsk+dwLNgA0SK99KrMyUWJknsbINGAZ8l1oCg7RHse+Eu/HsOuCzgE2HgXGqLM5K1CJSHngPuAhooKr7RORi3OCPjYGywFyguqompLUvq4mZ/tP7M/y74czrM48WlVr4HU7EGDkSbr3VJbIJEyA62u+ITDCFQk0skMRKJaAMcBw3/NSx9MqkNYr9VFVuBTqocmvAMjCrCczzKvAobt6YJJ2Acap6UlW34UbLbxyEY5kI9/LVL1OlWBX6TOnDoZOH/A4nYvTtC6+/DpMnwx132DiLJtspUFBjdCWg3uzOaUqrOTFpRLWeIryRfMlKlCJyHfCLqq5KtqkcsCPg+U5vnTFpKpinIB91+Yidh3Zy/+f3+x1ORBk4EJ55BkaNggcecMNVGRNsEit3AhNxTYgA5+N6QKYprcaBpJHq4zIVUNo3yf0LuDqlYimsS/FPRkTuAu4CyJMnT2ZCNBHmsvMv44krnuC5hc/R6cJOdL6os98hRYzBg938Y6+8AkWLuqRmTJANwLW8LQPQGN0ssVIqvUKpJjFVpnk/MzX3hWrKw4eISG3gAmCViIDLtitEpDGu5hU4/ML5wK+p7H84rkcLBQsWtP8NDQBPNX+KmZtncte0u2hSvgmlCqb7N2AyQAReesmN7PHss1CkCDz8sN9RmQhzUmP0lMS6uozESjSpVGICpdWxY1paO1DluszFmfw48hPQ0OvYURMYw18dO74EqlnHDnMm1u1ZR4PhDbim6jVM6TYF758lEwQJCdCjB3zyCQwfDnfe6XdEJitCqWOHxMoLwB9Ab+A+4B7gB43RJ9Mql1Zz4ktBiy6DVHWdiEzAzVkWDwxIL4EZk1zNUjX5T6v/8PCchxm5ciS31rPp74IlVy4YPRqOHIF+/SBvXujd2++oTIR4HLgdNyBxP2Amrgd7mjI0ALAIeXBd4RXYqEpITWpuNTGTXKIm0nJUS1bsWsHqu1dTqWglv0OKKMePQ8eObrDgjz+G7jaZQFgKpZpYIImVYsD5GqOr031tBkbsaA8Mxc3tIrjrWf1UmRWEWIPCkphJyfY/tlP7ndrUK1OPeb3nkSsql98hRZRjx9wULosWwfjxcP31fkdkzlQoJTGJlfm4+5KjgZXAXuBrjdGH0iqX0QGAr1KlhSrNgatw93gZE9IqFq3IG9e+wYLtC2yQ4GxQoABMnw6XXupqYp99ln4ZY9JQRGP0ENAV+EBjtAGQ7vxiGUlie1TZEvD8R9y00caEvD6X9KHzRZ3517x/sXbPWr/DiTiFCsHMmVC/vpshelbItM+YMBQtsVIGN0vK9IwWykgSWyfCTBH6itAHmAZ8mzSWYiaDNeasEBGGdRhG0XxF6TWpF6cSQupybkQoUgQ+/xxq1XLDU82d63dEJkw9A8wGtmiMfusN/rs5vUIZuSb2QRqbVZXbzijMbGDXxEx6Ptv4GZ3GdeJfV/yL51o953c4EWn/fmjZ0s0OPXOmm87FhLZQuiaWWRnqnRjqLImZjLh96u2MXDWShbcupEn5Jn6HE5H27IGrroLt213t7Ior/I7IpCWUkpjESj5cF/uaBMyUojGaZkUp3eZEES4Q4RURJonwWdKS5YiNOctebfsqFYpUoPfk3hw5dcTvcCJSqVLw5ZdQrpzrufjNN35HZMLIR7ihCq8BvsaN2HQ4vUIZuSY2BfgJeBPXUzFpMSasFM5bmFGdR/HjgR95ZM4jfocTsUqXdvePlSsHbdvCggV+R2TCRFWN0X8DRzVGRwHtgXSna89IEjuhyhuqfKXK10lLVqM1xg9XVryShy9/mKHfDWXWZutKl13KlnWJrHx5VyObP9/viEwYSJrV+Q+JlVpAEaBSeoUy0rGjJ1ANmAP8Ob+rKisyG2mw2TUxcyZOxJ+g0buN2H9sP2vuXkPxAsX9Dili/fab6+yxbRtMmwatWvkdkQkUYtfE7gA+BeoAHwDnAE9pjA5Ns1wGktjzwC24ETuSpsRTVVpmNehgsSRmztTK3Stp/G5j2ldvz6SbJtkgwdlozx5o3dr1Wpw6Fa5OaRIm44tQSmKZlZEktgGoE2rjJQayJGYy47Wlr/Hg7Ad5sc2LDGoyyO9wItq+fS6RbdgAU6a4a2XGf6GUxCRW8gLX45oQ/xycXmM0zdnrMnJNbBVQNAuxGROS7r/0fm64+AYen/s4C7Zb74PsVKKE67V48cXQqRPMmOF3RCYETQU64WYwORqwpCkjNbH5uDbKb/nrmpiq0ikLwQaV1cRMZh06eYhG7zbi0MlDfN/ve0qfk9Jk5CZYDhyANm1g9Wr49FM3Er7xT4jVxNZqjNY603IZqYnFAF2A/wCvAMuBqmd6IGNCUeG8hZl440QOnjhI94ndiU+M9zukiHbuuW5Yqnr1oGtXmDjR74hMCFkssZJul/rkMjqfWF2gJ25gxm3AJFXePNODZReriZmsGr16NLdMvoV7G93Lm+1C5qMdsQ4ehPbtYckSGDEC+vTxO6KcKSs1MRGpi5umKx+uCfAeVV1+xvuJlTW4uSqjcT3hf8S1+gmgGqN10iqf6szOIlQHugM9gP3AeEBUuepMgzQm1PWq04tVu1fx0pKXqFGyBvc0usfvkCJakSIwezZ07gx9+8LRo3CPveXh5gUgVlVniUg773mLTOynQ1aCSDWJARuAhUDHpKlYRHgwKwczJpQNaT2EDfs3MHDWQKoVq0abKm38DimiFSzo7h3r1g0GDIDDh+Gxx/yOypwBBQp7j4sAv2ZyP78B/XGXqdYA72uMZrhdP61rYtcDu4GvRHhXhFa46p0xESlXVC7GdB1DjZI1uPGTG9mwb4PfIUW8fPncdbEePeDxx2HwYIiAMcnDSbSIxAUsd51B2QeAF0VkB/AS8EQmYxgFNMQlsGs5w2ENM9I7sSDQGdes2NI74GRV5mQiWG+f8jRwJ276aYB/qepMb9sTuJGME4CBqjo7vf3ZNTETTD/98RON321MwTwFWXTbIsoWKut3SBEvIQH694f33oOBA+HVVyEqI93OTJakd01MRObiBuVN7kmgFfC1qn4qIjcBd6lqujMx/+MYsbJGY7S29zgaWK4xWj/D5c9kKhYRigE3At2yMmKHl8SOqOpLydZfDIwFGgNlgblAdVVNSGt/lsRMsMX9GsdVo66iUtFKLOi7gHPzn+t3SBFPFR56CF57DW67DYYPh1y5/I4qsmWxY8dBoKiqqrghbw6qauH0yv1jP7GyIjBpJX+enrSuif2DKr8Dw7wlO3QCxqnqSWCbiGzBJbQl2XQ8Y1LUsGxDpnSbQrsx7eg4tiNzbplDgdwF/A4roonAK69AoULw7LOus8dHH0Hu3H5HZlLxK9AcmI9rpUt3FuZUXCKxcsh7LEB+73lS78Q0E+MZJbEgu1dEegNxwMOqegAoBywNeM1Ob90/eG23dwHkyZMnm0M1OVGryq0Y3WU03SZ2o9vEbnx606fkyWWftewkAs884xLZo4/CkSMwYQIUsP8fQtGdwOsiEg2cwPs+PlMao1mqb2fbzM7ptKUuBfbherc8C5RR1dtE5H/AElUd7e3jfWCmqn6a1rGsOdFkp2Fxw+g/oz+dL+rM+BvGWyI7S4YNg7vvhiZNXC/Gc61FN+hCacSOzMq2mlhGL/CJyLvAdO/pTqB8wObzyXy3TWOCol/DfsQnxnPvrHvpNrGbJbKzpF8/KFYMbr4ZrrzS3VdW1vrYmGR86f8jImUCnnYB1nqPPwO6i0heEbkAd/f2Gd8BbkywDWg8gLeufYspG6bQbWI3TiWE7KQOEeXGG2HmTPjpJ1cj27TJ74hMqPGrE+sLIrJGRFYDV4G7iVpV1wETgB+Az4EB6fVMNOZsCUxkncZ14ugpa8I+G1q3drNEHz0KV1wB333nd0QmlGTbNbGzya6JmbPpvRXv0W96PxqVbcSMnjNsZuizZNMmN6Hm77+7Oclahsy0vOErEq6J2e2ExpyhO+rfwac3fcrK3Stp9kEzdhzc4XdIOUL16rBoEVSoANde66ZyMcaSmDGZ0PmizszuNZtfDv9CkxFNWLl7pd8h5QjlysGCBdCwobteNiy77lg1YcOSmDGZ1LxScxb0dTNCNx3RlMnrJ/scUc5QrBh88YWrjfXvD7GxNt5iTmZJzJgsuKT0JSy/Yzm1StWi64Su/Gfhf4iE68yhrkABd12sd294+mm48044fdrvqIwfLIkZk0VlCpVhfp/59KjVgyfnPcnNk262notnQe7cMHKkG/n+/ffhuuvcCB8mZ7HeicYEiary/DfPM3jeYGqUrMEnN37CxSUv9jusHOHdd93oHpdcAjNmQOmUxgoy/xAJvRMtiRkTZF/++CU9J/XkyKkjDO8wnJvr3Ox3SDnCzJmus0epUjBrFlx0kd8Rhb5ISGLWnGhMkLWq3Irv+31Pw7IN6TW5F/2m9bPmxbOgXTv4+ms4dsyN7vHNN35HZM4GS2LGZIOyhcryZe8vebzp47y74l3qDavHsp3L/A4r4jVsCEuWQMmSbqSPiRP9jshkN0tixmST6Khonm/9PPP6zONkwkmajmhKzFcxnE6wbnTZqXJlWLwYGjSAm25yk2yayGXXxIw5Cw6eOMjAzwfy4aoPaVi2IR90+oBapWr5HVZEO34cbrnFjewxcKCbcNNmiv47uyZmjMmQIvmKMKrzKCbeOJGf/viJesPqMXjeYE7En/A7tIiVPz+MHw8PPghvvAGdOsHhw35HZYLNamLGnGX7ju3j4TkP8+GqD6levDrDOwyneaXmfocV0YYOhXvvhZo1Yfp0KF8+/TI5gdXEjDFnrESBEozqPIo5veYQnxhPi1EtuH3q7ew5usfv0CJW//5/zUvWuDHExfkdkQkWS2LG+KRNlTasuXsNjzZ5lA9Xf0j1N6vz6pJXbcLNbHL11a7DR758bqboSZP8jsgEgzUnGhMCNu7byIOzH2TWlllcWPxCXmv7Gm2rtvU7rIi0Z4+7PrZ0KQwZAo8+CiJ+R+UPa040xgTFhSUuZObNM5neYzqJmsi1H19L+zHtWf3bar9DizilSsG8edCtGzz+ONxxB5yyym/YsiRmTAhpX709a+9Zy4ttXmTxjsXUHVqXWybfwo8HfvQ7tIiSPz+MGQP//jeMGAFt27oZo034seZEY0LUgeMHeGHRC7y+7HXiE+Pp16Afg68czHnnnOd3aBFl9Gi4/XaoWBE++yxnjblozYlZICL3ichGEVknIi8ErH9CRLZ4267xKz5j/HZu/nN5vvXzbBm4hdvq3cY7ce9Q+Y3KDJoziN1HdvsdXsTo1Qu+/BL++AMuvdQNHmzChy81MRG5CngSaK+qJ0WklKruEZGLgbFAY6AsMBeorqoJae3PamImJ9i8fzPPLHiGMWvGkCdXHu6sfyePNn2U8wuf73doEeHnn12Hj1Wr4IUX4OGHI7/Dh9XEMu9uYIiqngRQ1aQbZDoB41T1pKpuA7bgEpoxOV614tX4qMtHbLx3Iz1r9XQ1s9cr029aP7Yd2OZ3eGGvQgU38v0NN8Ajj7hZo0/YgCqpEpEbvZa0RBFpmGzbWWtR8yuJVQeaicgyEflaRBp568sBOwJet9Nb9w8icpeIxIlIXHx8fDaHa0zoqFqsKu93ep/N923m9nq3M3LVSKq+WZUbP7mRJTuW+B1eWCtY0A1V9eyz7lpZ8+bw669+RxWy1gJdgQWBK70Wte5ATaAt8LaIZNuoldmWxERkroisTWHpBEQD5wKXAY8AE0REgJQq7ym2d6rqcFVtqKoNo6Ojs+s0jAlZlYpW4p0O7/DjwB95pMkjzP1xLk1GNKHJ+02Y+MNEEhLTbIU3qRCBwYPdzdDr1rnpXZYv9zuq0KOq61V1YwqbzmqLWrYlMVVtraq1Ulim4mpYk9RZDiQCJbz1gaOanQ/Y/0HGpKFc4XIMaT2EHQ/u4I22b/Db0d+48ZMbqfZmNV5d8ioHjh/wO8Sw1KWLm5ssb143wsdHH/kdUbaITmrR8pa7grDPDLeoBYNfzYlTgJYAIlIdyAPsAz4DuotIXhG5AKgG2P9AxmTAOXnO4b5L72PTvZuYdNMkyhYqy0NzHqLcK+W4beptLP9lOZFwS83ZVLs2fPstXH65u0b26KOQEFkV3PikFi1vGR64MZ0WtdRkuEUtGPxqhxsBjBCRtcApoI+6v651IjIB+AGIBwak1zPRGPN3uaJy0aVGF7rU6MLK3SsZGjeU0atH88HKD6hfpj53N7ybHrV6UDBPWHdKO2tKlIA5c+CBB+DFF13vxTFjoHhxvyPLfqraOhPFzmqLmt3sbEwOcOjkIUavHs07ce+wds9aCuctTI9aPehbty+XlrsUifS+5EHy3nswYACULeuumdWr53dEWROMLvYiMh8YpKpx3vOawBj+ulXqS6BadlVILIkZk4OoKot3LGbYd8OY+MNEjscf58LiF9Lnkj7ccsktds9ZBixfDtdfD/v2wbvvupulw1VWkpiIdAHeBEoCfwArVfUab9uTwG24FrUHVDXbbiG3JGZMDnXo5CEm/jCRUatGsWD7AgShTZU29LmkD50u7GTNjWnYswduugm+/hruuw9efhly5/Y7qjMXCTc7WxIzxrD19618uOpDRq0axfaD2ymQuwAdq3eke63utK3alnzR+fwOMeTEx8Njj8Err8AVV8Ann0Dp0n5HdWYsiYUIS2LGBEeiJrJw+0LGrR3HxPUT2XdsH4XzFqbzRZ3pXrM7rSu3JneuMKxyZKNx49wAwkWKwKefup6M4cKSWIiwJGZM8MUnxjNv2zzGrR3HpPWTOHjyIMXyF6PLRV3oclEXWlVuZTU0z+rV0LWrG3/x9dehf//wGHfRkliIsCRmTPY6GX+SOVvnMG7dOKZtnMbhU4cpmLsg11a7ls4XdqZ99fYUzVfU7zB9deCA6+Qxcyb07Qtvv+3mLQtllsRChCUxY86ek/Enmf/TfCZvmMzUjVPZfWQ30VHRXFXpKjpf1JmO1TtSvkj59HcUgRIT4ZlnIDYW6tZ118mqVvU7qtRZEgsRlsSM8UeiJrL8l+VM2TCFyRsms2n/JgBqlarFtVWv5dqq19K0QlPy5Mrjc6Rn18yZcMstrvPHBx+4psZQZEksRFgSMyY0rN+7nhmbZzBryywWbl/I6cTTFMpTiFaVW/2Z1HJKLW37dtcNf/lyN9rHf/8LeUIsl1sSCxGWxIwJPYdPHmbetnnM2jKLWVtm8fPBnwGoWbImbSq3oVXlVlxZ8UoK5y3sc6TZ59QpGDQI3nzT9VocPx7Kh1AOtyQWIiyJGRPaVJX1+9Yza/MsPt/6Od/8/A0n4k+QS3LRqFwjWlZqSavKrWhSvklE9nj85BPXDT9PHjdPWdu2fkfkWBILEZbEjAkvJ+JPsGTHEuZtm8eX275k+S/LSdAE8ubKS9MKTWlZqSUtKrWgYdmG5I3O63e4QbFpk5s1eu1aN19ZTAzkyrapIjPGkliIsCRmTHg7dPIQC7cv5MttXzJv2zxW/bYKgLy58tK4XGOaVWhGs4rNuPz8yymSr4jP0WbesWNw772us0erVvDxx3Deef7FY0ksRFgSMyay7D26l29+/oZvfv6GhT8vZMWuFSRoAlESRZ3z6rikVqEZV1S4gjKFyvgd7hkbMcKNhn/uuTB2LDRv7k8clsRChCUxYyLb0VNHWbpzKQt/Xsg3P3/Dkp1LOHb6GACVilbi0nKXuuX8S6lXuh75c4f4Xca4UT5uuAG2boWnnnJNjGe7edGSWIiwJGZMznI64TQrd69k4c8LWbpzKct+WfZn78foqGguOe+SP5PapeUupVrxakSJXxPZp+7wYbjnHtfZo3lz17xYrtzZO74lsRBhScwYs/vIbpbtXMayX9zy7S/fcvjUYQCK5itK43KNaVCmAfXL1Kd+mfpcUPSCkJkM9MMPXTLLlw9GjYL27c/OcS2JhQhLYsaY5BISE9iwb4NLajuXsfzX5azds5b4xHgAiuQtQr0y9ahfuv6fia168erkivKny+DGjdCtG6xaBQ8+CEOGZP/N0ZbEMntQkfHAhd7TosAfqlrX2/YEcDuQAAxU1dnp7c+SmDEmI07Gn2TtnrWs2LXCLbtXsGr3Kk4mnASgQO4C1C1dl/ql61O3dF1qn1ebmiVrnrUJQk+cgEcegbfeggYN3DQv2Tn2oiWxYAQg8jJwUFWfEZGLgbFAY6AsMBeorqoJae3DkpgxJrNOJ5xmw74Nfya273d/z/e7v+fIqSMACELlcytT57w61C5Vm9rn1aZ2qdpULVY122ptU6bAbbe5sReHDYMePbLlMJbEsnxw1yD9M9BSVTd7tTBU9Xlv+2zgaVVdktZ+LIkZY4IpURP58cCPrPltDWv2eMtva9j8+2YSNRGAfNH5qFmy5p9JrXap2tQoWYNyhcoF5Vrbzz9Dz56waJFLaG+8AQWDnG4siWX14CJXAq+oakPv+VvAUlUd7T1/H5ilqhNTKHsXcBdAnjx5Gpw8efLsBW6MyZGOnz7OD3t/+DOpJSW43Ud2//maQnkKcVGJi6hRsgY1SrjlohIXUaVYFaKjos/oePHxblqX556DCy9095TVrRu887EkltaOReYCpVPY9KSqTvVe8w6wRVVf9p7/D1iSLInNVNVP0zqW1cSMMX7ae3Qva/esZf2+9WzYt4H1+9azfu96fjn8y5+vyR2Vm2rFq/2Z1GqUqEGNkjWoXrw65+Q5J839z5vnpnbZtw+ef96Nih8VhDsGLIll5cAi0cAvQANV3emts+ZEY0zEOHTyEBv2bXCJbe96l9z2rWfr71tJCLjUX/qc0lQrVs0txatRtVhVqhVzP5M6lezfD3fc4a6XXX01jBwJZbI4WIklsawcWKQt8ISqNg9YVxMYw18dO74EqlnHDmNMJDmVcIotv29h/d71bNq/ic2/b2bL71vY/PvmvzVNApQtVDYgqVXjp++qMvKVapxzqiofDC9Ax46Zj8OSWFYOLDISd/1raLL1TwK3AfHAA6o6K719WRIzxkSKwycP/5nQNu/fzJYDW9i8fzObf9/MnqN7kr24NA3y9CTuPy9n6liWxEKEJTFjTE5w8MRBtvy+xdXi9mzm03k/ctG5dfjk4Qcytb+sJDERuRF4GqgBNFbVOG99G2AIkAc4BTyiqvMyFWBG4rAkZowxOVMWk1gNIBEYBgwKSGL1gN9U9VcRqQXMVtVsGxHyzPp7GmOMMYCqrgf+cU+cqn4f8HQdkE9E8qpqttwHZUnMGGNyrmgRiQt4PlxVhwdx/9cD32dXAgNLYsYYk5PFJw02kZKM3O+bRtmawH+Bq7MWYtosiRljjEmRqrbOTDkROR+YDPRW1a3BjervQm+WOGOMMWFLRIoCM3D3AS/K7uNZEjPGGHPGRKSLiOwELgdmeCMsAdwLVAX+LSIrvaVUtsVhXeyNMSZnioSbna0mZowxJmxFRE1MRBKB41nYRTRumKucIqedL9g55xR2zmcmv6qGdWUmIpJYVolIXFrdTCNNTjtfsHPOKeycc56wzsDGGGNyNktixhhjwpYlMSeYw6yEg5x2vmDnnFPYOecwdk3MGGNM2LKamDHGmLBlScwYY0zYytFJTETaishGEdkiIo/7HU+wiEh5EflKRNaLyDoRud9bX0xEvhCRzd7PcwPKPOG9DxtF5Br/os88EcklIt+LyHTveaSfb1ERmSgiG7zf9eU54Jwf9D7Ta0VkrIjki7RzFpERIrJHRNYGrDvjcxSRBiKyxtv2hiSf+CtSqGqOXIBcwFagMm4a7VXAxX7HFaRzKwPU9x4XAjYBFwMvAI976x8H/us9vtg7/7zABd77ksvv88jEeT8EjAGme88j/XxHAXd4j/MARSP5nIFywDbcDboAE4C+kXbOwJVAfWBtwLozPkdgOW5cQwFmAdf6fW7ZseTkmlhjYIuq/qiqp4BxQCefYwoKVd2lqiu8x4eB9bgvgE64Lz68n529x52Acap6UlW3AVtw70/Y8KZ+aA+8F7A6ks+3MO7L7n0AVT2lqn8QwefsiQbyi0g0UAD4lQg7Z1VdAPyebPUZnaOIlAEKq+oSdRntw4AyESUnJ7FywI6A5zu9dRFFRCoB9YBlwHmqugtcogOSRpaOhPfiNeBRIDFgXSSfb2VgL/CB14T6nogUJILPWVV/AV4CfgZ2AQdVdQ4RfM4BzvQcy3mPk6+PODk5iaXUPhxR9xuIyDnAp8ADqnoorZemsC5s3gsR6QDsUdXvMlokhXVhc76eaFyT0zuqWg84imtmSk3Yn7N3HagTrtmsLFBQRHqlVSSFdWF1zhmQ2jnmhHMHcnYS2wmUD3h+Pq5pIiKISG5cAvtYVSd5q3/zmhnwfu7x1of7e9EUuE5EfsI1C7cUkdFE7vmCO4edqrrMez4Rl9Qi+ZxbA9tUda+qngYmAU2I7HNOcqbnuNN7nHx9xMnJSexboJqIXCAieYDuwGc+xxQUXi+k94H1qvpKwKbPgD7e4z7A1ID13UUkr4hcAFTDXRQOC6r6hKqer6qVcL/Hearaiwg9XwBV3Q3sEJELvVWtgB+I4HPGNSNeJiIFvM94K9z13kg+5yRndI5ek+NhEbnMe696B5SJLH73LPFzAdrheu5tBZ70O54gntcVuKaD1cBKb2kHFAe+BDZ7P4sFlHnSex82Esa9mIAW/NU7MaLPF6gLxHm/5ynAuTngnGOBDcBa4CNcr7yIOmdgLO6a32lcjer2zJwj0NB7n7YCb+GN0BRpiw07ZYwxJmzl5OZEY4wxYc6SmDHGmLBlScwYY0zYsiRmjDEmbFkSM8YYE7YsiZmIJCIJIrIyYKnkd0xZISKdReQp7/HTIjIo2fafRKREGuXnBo58bkykiPY7AGOyyXFVrZvSBu/mT1HVxJS2h6hHgeuyUP4j4B7gueCEY0xosJqYyRFEpJI359bbwAqgvIg8IiLfishqEYkNeO2T3txMc705qwZ56+eLSEPvcQlvmKukecxeDNhXP299C69M0pxfHyfN6SQijURksYisEpHlIlJIRBaKSN2AOBaJSB0RqQ6cVNV9GTjP/gG1z20i8pW36TOgRxDeSmNCiiUxE6nyB3yZT/bWXQh8qG7A3AtxQ/Q0xo180UBErhSRBrihq+oBXYFGGTjW7bgR1Rt5r7/TGwIIbz8P4OZ9qgw09YY5Gw/cr6qX4MYEPI6bRqYvgJe48qrqatzYkCuSHfPBwOZS3IC4qOpQrwbaCDfawyve+gNAXhEpnoHzMSZsWHOiiVR/a070roltV9Wl3qqrveV77/k5uKRWCJisqse8chkZT/NqoI6I3OA9L+Lt6xRuHLud3r5WApWAg8AuVf0WQL0ZBkTkE+DfIvIIcBsw0ttfGdy0K4FeVdWXAs7vp2TbX8eNITktYN0eXLLbn4FzMiYsWBIzOcnRgMcCPK+qwwJfICIPkPqUFfH81XqRL9m+7lPV2cn21QI4GbAqAfc3JykdQ1WPicgXuOlGbsKNfQeullYklZj+QUT6AhWBe5Ntyufty5iIYc2JJqeaDdzmzbmGiJQTkVLAAqCLiOQXkUJAx4AyPwENvMc3JNvX3d70N4hIdXETVKZmA1BWRBp5ry8kbqZicE2KbwDfqmrS7L7rgaoZOSmvOXQQ0Cuw44p3La60dw7GRAyriZkcSVXniEgNYInX1+II7ot/hYiMx438vx1YGFDsJWCCiNwCzAtY/x6umXCFlyz2ksZU8Kp6SkS6AW+KSH5c7ag1cERVvxORQ8AHAUUWAC+LiGj6I3bfCxQDvvLOK05V78Al36WqGp9OeWPCio1ib0waRORpXHJ5Kb3XBul4ZYH5wEXJalKvA9NUdW4m9/s68JmqfhmUQI0JEdacaEyIEJHewDLc3HbJ72H7D1AgC7tfawnMRCKriRljjAlbVhMzxhgTtiyJGWOMCVuWxIwxxoQtS2LGGGPCliUxY4wxYev/Adwurdeofw2NAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax1 = plt.subplots()\n",
    "ax1.set_title('Magnitude response(dB) and Phase Response')\n",
    "ax1.plot(w, 20 * np.log10(abs(h)), 'b')\n",
    "ax1.set_ylabel('Amplitude(dB)', color='b')\n",
    "ax1.set_xlabel('Frequency(Hz)')\n",
    "ax2 = ax1.twinx()\n",
    "angles = np.unwrap(np.angle(h))\n",
    "ax2.plot(w, angles, 'g')\n",
    "ax2.set_ylabel('Phase(radians)', color='g')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cc6eef62",
   "metadata": {},
   "source": [
    "#### 您也可以点击工具栏上的幅值和相位响应按钮，或选择分析 > 幅值和相位响应显示绘图。\n",
    "\n",
    "#### unwrap 函数在频率分析中也很有用。unwrap 根据需要对相位增减若干个 360° 以将其展开，使之在 360° 相位不连续点处保持连续。要了解 unwrap 的作用，请设计一个 25 阶低通 FIR 滤波器："
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 273,
   "id": "f05d9bca",
   "metadata": {},
   "outputs": [],
   "source": [
    "h = signal.firwin(25, 0.4)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "84db5725",
   "metadata": {},
   "source": [
    "#### 用 freqz 获得频率响应，并以度为单位绘制相位："
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 274,
   "id": "cc65f690",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, 'Phase(radians)')"
      ]
     },
     "execution_count": 274,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD4CAYAAAAQP7oXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABRVUlEQVR4nO29e5hcV3Un+lv9fneVpFbr0d2qkrBNbOIYrJiHEy4JdmwGEgO53Jg7IQzJjBMGMhkeMwMzN9NUZpwvHySeXL5cSAwmmBAwHsDYgHECJLEhMdjyA9tCGMuSWmo925b6oW5J3epe94+qXTp1ap9z9j5n71OnSvunr79uVZ/Ta5+qvfd67t8iZoaDg4ODg4MK2ho9AAcHBweH5oFTGg4ODg4OynBKw8HBwcFBGU5pODg4ODgowykNBwcHBwdldDR6ALaxYcMGLhQKjR6Gg4ODQ1Phsccee4GZR/yvt7zSKBQK2LVrV6OH4eDg4NBUIKIp2esuPOXg4ODgoAynNBwcHBwclOGUhoODg4ODMpzScHBwcHBQhlMaDg4ODg7KcErDwcHBwUEZTmk4ODg4OCjDKY0mwxNHn8DDhx5u9DAcHBwuUrT84b5Ww3//p/+OqdkpPPXupxo9FAcHh4sQztNoMqysrmDfqX2w2TzrzMoZvPyvXo7vTX3PmgwHBwd1/HD6h7jqL6/CwrmFRg/FKQ2TeO7F57BnZo9VGQzG4soiXlh6wZqMk2dO4sljT+K7+79rTYaDQxp4aOohnDpzqtHDSIzdM7vxo+M/wu6Z3Y0eilMaJvH+v38//s29/yYVWftn91v724yyF3Ng9oA1GQ4OtrHGa7juc9fhz3/w540eSmKIyEIW1qRTGgZxZuUM9p7ca1WGmDz7T1lUGkKGRcXk4GAba7yGlbUV7D1ld00+cfQJvOFv34Cz589akyEMOZvrXhVOaRgEg3HyzEnMn5u3KgNIx9PIwgR1cIiLNAwsAHjk8CN4YO8DVkPTWTLknNKwgFS8gBRkTM9PY3l12ZocB4c0YHujTdWQc0qjtZCGNZDm5GEwDs4dtC7HwcEGxFo5dvoYzqycsScnxXxDFrx/pzQMIs2wThqKCcjGJHVIH9Pz0zi/dr7Rw0gEb1m6zQ09jXUvnuXg3EGsrq1ak6MCpzQMIhVPoyJjanbK2uTxLrYsuMMO6eLc+XN46V+8FLc/dnujh5IINcZPCmsyDUNuZW0FhxcOW5OjAqc0DCLN2ObK2gqOLByxKgNwnkYr4IfTP8RH//mjytcvry5jcWURTx1vbtaBGuOn2b3/lJ8lDE5pWEAabipgb5I6T6O18KXdX8KHvvMhnDt/Tun6NV4D0FqffVrnmmwxNaTlNanAKQ2D8LqpNifP1sGtAFJKvLXQxnGxYo3XtIoaqkqjyb3MtMNTSytLOLF4wpocgUZ/Lk5pGISYpEsrS5hZmrEmZ2J4AgSyNnnEc/R29DZ8gjokh67nID7/qbmp6r3NiLRCOmkoJ/EsvR29DTfkGqo0iOgzRHSCiJ7xvPYRIjpMRE9Wvv6V53cfJqK9RPQsEd3QmFEHI41Jyszo7ujG1qGt1idoMV/EzNIMTi+ftiLHIR3onu0RimJ5dRlHF45aG5dtNCKkY9uQK+aLF7fSAPBZADdKXv9fzHxV5et+ACCiywHcDOCKyj2fIKL21EaqAAZjtH8UgL3QEYNBIBRz9iaPmKDb89sBlCu1HJoXQgmozslWyWmJ5xjtH8Xs2VnMnp21Kgew72lsz29vOP9UQ5UGMz8E4KTi5TcBuIuZzzHzfgB7AVxjbXAxUcwXAdidPERUtjgsh46KObvP4hAPx04fw92771a+XreqzxuSaoXwpDB+bBpyANDV3mXf08gVcXj+sHJRgw002tMIwnuJ6KlK+CpfeW0rgEOea6Yrr9WBiG4hol1EtGtmxl5uwQ9mxkDXAEb6RqxOHgKhMFywRvNRDU8JpdECG0cr4fNPfR6/8eXfUKbH181p1CiNJjYY/B6zzZAxAKvev0AxV2w4U0MWlcYnAewAcBWAowD+rPI6Sa6Vligx8+3MvJOZd46MjFgZpFSuCB1ZjjsKT8PW5BGLbXRgNBOJN4darKyuANDPUahen6XyziTwhnQAi94/7IeO0noWFWROaTDzcWZeZeY1AJ/ChRDUNIBxz6VjAOycbouJaujIZr4hBS9AyGijNhRyhabeOFoRuuEm8XmqFjW0SnhKvE/5njyGu4dTCefaovlIy2tSQeaUBhFt9vz3LQBEZdV9AG4mom4iKgK4BMAjaY8vDN4ktS2aD683A9ixOMQErXpNTbxxtCJ0PQevElCxhFstEU5EVo0frxdgi+ZDyNg6tBWdbZ0Xr6dBRF8E8DCAy4homoh+B8BHiehpInoKwC8BeB8AMPNuAHcD+DGABwC8h5kby9wVgEKuYI3mQ3gzWwcrk8fihu71mmz2JHfQg3aOAnqeg/j7Wwe3Ynp+uhoOa2bYDBmn4QUIGe3Ujm25bRev0mDmtzPzZmbuZOYxZr6Dmd/BzD/LzFcy868x81HP9bcy8w5mvoyZv9XIscvgrWwC7HkBBEJ7WzsmhifsyPAoiGKuiPlz8zh19pRxOQ7xoEvFzczo7egFoDYnhdLYsW4H1nitaenxazzmXNEazYf3XBNg3zsr5hrr/WcuPNXM8IanAPtxR1vWkz88BTR3bLvVEKcaamP/RvR19il9jnWWc5OGqLzhqWKuaI3mQ7xf24a3WWNq8D/LRetptBqEp1Gl+bDkBRCVC8mE9WRDBnBhggLNu3G0IryH9VRoPhiM9rZ25bh+1dPI7wDQvAaD1PixOI9tMjX4n+WFpRcaxtTglIZhEMj65CFcUBonFk9gcXnRuAzAeRpZhfh8VGk+1ngNbdSmbKEKpTE+NI52am/4CeSkqDF+LHoBgL2zGlJDrkFr0ikNg/DWt9uKO9Z4GpUN3daiJiLkenLI9eSafuNoJegevlvjNa24vvh9V3uXtbxZGvA+ZyFXAND81YZeQ65Ra9IpDYPwb+hWQkc+TwMwvxD8m0qjY6gOtdAlxmTmsqeRVytqEEqpSlfTpJ+9dzPv7+rHxv6NqYRzjywcMU7z4TdIgcaFjJ3SMAj/hm6L5kOgaj0Ztmy8iw3IBrOmwwWs8RraK1ydyp6GRlhDfP7VkFaThia9mzlg3/gR3hyDMTVnluTT+ywb+jYoFzXYgFMaBuFPUtug+fDKqFbEWPI0/Al3d1YjG1jjNXS1d2HL4Balz55xwdMAohWN8DQEI8DxxeNYWllKPvCUITV+LJ6hEDIA86Ej77M0uoLKKQ3D8E5QwI4XIGTYOunqX2yFXAFnz5/FsdPHjMpxiIeqElD0AryJcCB6TlbDU57ycZ1N8KnjT+EzT3xG+Xrb8Bo/Nmg+0khS1xlyDfT+ndIwiDTijl5PQ8hxFOkXF6rhJsWNQyTCh3uGke/JR97j5R6LY/x89snP4pav34Lza+eV77EBv2csmBpM03x4jawtg1us0nx4w9/7TzWGqcEpDYPwbujVyWNhQxcTB7gQpzU5eWRWDeDKbrMCr+egQvMhEuGAmoVakwiPYTCsrq1ilVcxPT+tfI8N1IWnUvAC2tvs0Hx4DVKg/CwLyws4eUa1HZE5OKVhEN7QkS2aj7rJo1gRE0eGNzwFOE8jK6hWQ+WKSjQfQskA5c9SJxG+aWATejp6tDZaXUJFWwg0flLyAkwi7WcJg1MaBlEXOrIQd5SFpwCzC9Q/Qfs6+zDaP9rwTcChjOq5C8WkqwhnAReKGsJOknsT4XHyZrrU7bbgN36qTA2Wqg0FbCSp0/KaVOCUhmHUhY4sJsIBuwf8/HIavQk4lOFNhAPRm7O4HijPyXOr50KLGryJcHGPzmefFU9DQCjMrvYujA2NWa82LOQKxmk+nKfRopBZHKqNb5RlBHkaBieP/zmEHKc07OLQ3CE8fOjhyOuE5zA2NIaOtg6laiidqj5vIhzQ5zgT9zd6vsjyfFa8f0lpL2BHaQoZQ91DWNe7znkazQ5ZeAow7wV4PYBqRYzF8BRQ3jgOzR1qeEVMK+NP/+VP8cYvvDGyqEHkKFTzZjWJcIUSWm8iHCjP49mzs5g9O6v0HLosvLbg38yBlPINKRpyB+YOGJOhCqc0DKIudBSjxl1Fhh+mrSfpYssXM1ER08o4e/4sTp09FVnU4FcCKtVQ3kQ4EL6heRPhQgZgrye5LQQZPzZoPryw4WlIn6VBXTWd0jCIIE/D9OTxygBg/IBf0GIDGr8RtDKqCWSNcJNKNZQ3Ed7b2YtNA5tC7/EmwgH9+Ll4jqOnj+LMyhmle2wgyPgxzdTgN+RG+kaMMzXInqUwXFCmxzcJpzQMw/uh2po8XhmAPZoPWcK90SGHVoZqWMfrORRzxUiaD28iXNwTJkOWCAfi9SQ3zcEUB9ZzgFy7JsX5FqMRhgBPI6qowQYa3SP8M0R0goie8by2joi+TUTPVb7nPb/7MBHtJaJnieiGxow6GH6LwwZHjMzTKOaKRmk+ZCGw8aFxtFGb8zQsQjWsU1MNpZA383om4p7Q8JQvEZ7ryWGoe0jb04gal20EJcIBw94/JGsyhWrDRnn/jfY0PgvgRt9rHwLwXWa+BMB3K/8HEV0O4GYAV1Tu+QRRheozI5Bu6BbijnWehmEvQGbVdLZ3WilXdLgA1fMN/nMXQHQ1lN/TCCtq8IendI2fNV5Dd3t35LhsQxbS2Tyw2TjNh9/TAMzTfATlMoH0vf+GKg1mfgiA/xz8TQDurPx8J4A3e16/i5nPMfN+AHsBXJPGOFUhCx0VhgtGaT6CqigAcwtUttiEHKc07EE1PCXzNKLCTX6lscqrODR3KHQccZOua7yGzYOb0d3e3dD5IjN+bNF8yLx/kzQfMsVkqzVCFBrtacgwysxHAaDyfWPl9a0AvLN8uvJaHYjoFiLaRUS7ZmZmrA7WiyBPwyTNR1AiHLDraQCNq9a4WFA936CRCB/tH0VvR29kYltaoBEwX/zVU4Be3oyZ0U7qPcltIdT4MRye8sO49y8JgfV09GDzwOaLy9PQBElek85gZr6dmXcy886RkRHLwwqHDS/AvwhUKmLiQLbYGl0R08oQFr4KzYc3dBS1OfsT4VEWqj8RDpQ/+zPnz+D44nGl5xC9OLJoZFjJM0rWCmBw3UtkAI1hasii0jhORJsBoPL9ROX1aQDjnuvGABxJeWyhkFY2Wcg3+C0OwOxCkFlOQgaQjYqYVoTYrKMqYupyFBGthf2J8GpRQ5CnwRJPQyOB7O0UmLXwFFB+FpM0HzIvIC2Sz0Z0Vsyi0rgPwDsrP78TwL2e128mom4iKgK4BMAjDRhfIIIqmwCzcccgi8NUpUrYYgPcWQ1b8CprrXCTQgmtVwF0tndifGg88B5/IlzIANQ2QW+nwJNnTmL+3HzkPQIPTT2ELz3zJeXro8YByD1mwO6aNM3UIFNMQKWoYf5QJD2+STS65PaLAB4GcBkRTRPR7wD4EwDXE9FzAK6v/B/MvBvA3QB+DOABAO9hZrMtuBJC5mmoNr7RkSGD6EpmguYjcrG5ZLgVeHt/R3kO/g09jObD75kA4UaGLBEuLGcVw6TKwhtjc/6LR/4CH/j7DyhfH4ZI48dwDtAPk6GjsPDUGq+lytTQ6OqptzPzZmbuZOYxZr6DmV9k5tcz8yWV7yc919/KzDuY+TJm/lYjxy5DYOjIsBcQZHGYpvnwyxEVMY2svW9lrPEaJoYnAETTfOiEjvyeCRAe1pAlwvu7+rGxf6OSAqj2+4ixOa/yqjGaj7Q8jTAvwGYiHGhMr5sshqdaDsYnj8TiMFl+F2Q5tVGblXJFhzKYGX2dfeWKmKjwlITjLKwaqs7TCClqkCXCxT2q7WV1epL77zVN8+HHhr4N6O/st+8FKPQuMSEDSDdk7JSGQQRt6CZpPsK8GcCMxRFkoQGNSbxdLFDt/e0PN6lUQ8nCU4C8qEGWCBf3KPckJ8K63nUY7BrUa+BkkFY9KDwl3mNj3n+Qp5EvYnl12SrNx/jwONqp3XkazYqwDd0kzYcMJmk+ghYb4A742YS3uVJUYtv72eR78xjuHg5NbOuEaGSJcHHPwbmDWF0LTyWK51BRgEGyjczjKOPH4Dy27QUEKaaOtg6MDwcXNdiAUxoGEeZpAOa8AJmMqIoYXRmAfCEUcgXtihgHNYjNvZArRNJ86HgBQYlwQD4nZYlwoPzZn187H5k38xMq6myaJlvFhhk/4gyJKe9fBqPef0B4Ckjf+3dKwzDSyDfIFgFg/qBPaBjMhaiMw5sLCKP5CMpR6CTCNw1sKtN8SO6RJcKFDECB5sSzwemGZm00cAoy5EzRfEQmqS16GkD63r9TGgYRtDBMVjgEeRqAOYsjbIHH8ZoePPAgrvvcdanWkjcjVKuOgsJNQZuzTMlUT2yHeRpBB1UVaE683FiLK4t4YemF0Hv8sk2Gp2RIwwswTfMRuO7zRRw7fSw1pganNAwiyBowSfMR5aaaoPkIjQXH8DQePfIovrv/uy4XEgHvSWpAP7EdRPMhu17cE6Y0/PdMDE+AQFosvLoGUxqJcCCdfANgzvtXMeTSKoV3SsMgIuOOtikFQipidBC22Nb3rsdA10C8ihgX0gqF2NyjKmKkOYqIxLZOJVxQ9VRXe5cSPX4NC2/MBk4maD6UjB/bXoAp7z9CMQHpndVwSsMgUrE4ImKbQPLNOWyxxeEUshGnbkWIzVZUxISd2A465Sy7R6ZkgLIXcOrsKcydnav7+0BwTks3PAXot4oFDMzjEONnqHsI63rXpRLONUHzEWWQAukZZU5pGEbYBxtWEaOK0Mmj0MVNB6EKUPPAFuA8jSh4PYIwxSzLUYSFgWRKRsgA6udLUCI8alzV+z1zdKBrABv6Nmh5GtUGTml4AZYNuUKugDVew6F5eVGDCRmbBjahp6PHeRrNiCiLwwTNR9TkMdH4JqrSRbcixmQZZSvD6xFEVUP5N/S+zj6M9o8GVkMF5TSA+s8lKKchxhVF8yHjxtLp+rc9v708LkMecxBM5huiDDkTBlOQDCLCtuFtLqfRjFCKO1qcPKZoPsLCU0B5E4hTEeM4q8JR08Y1H07zodNbITARHhDWCKqeEjKA8LxZndLQ2JyZGSP9I0ZoPsLCU4BZmo/IkLFtQy7FvhpOaRiEUtwxBS/AZiwYiBGnNlgR08rwH4oDgmk+dEJHQUomiOYjKBHuHVdoT3Kf8VTMFTE1OxV5klyMtZ3ajdB8qBg/Jmg+wjyaalGDAa8paD0C6R7wc0rDIMI+2FQnTwqxYKAxFTGtDGnv76BqqIDDZDKajyAlE0TzEZUIB/R7kq+sreDIQnTPNJMNnJSNHwNGVtBaqRY1zB2wJgMov8eyogYb0FYaVKJ+KlVI/x3qEDV5TCyE0MkTo/FNnYyIWLBu7b3X/XfJ8GB4PYKoxHaQFyCj+Qi6Xtzj/0zCEuFbBregs61Tj4VXwzOt4d8yRPMRtF5MHbpNwwuIlJFi2W2k0qAStVGJ/m8q0TepRCcA/ATAUSrRbirRx6hEl1gfZZMgaoJH9XJWkqEwQYFkm3OUhTbYPYj1veuVZdSUUV7EIaqo2Ll3c4+i+dBNbEdVEHnnblgiXCVvVteOVmNOVhs45ZPTfCgbPxY9DcCc9x8lA0jHKFPxNP4RwA4AHwawiSd5nCd5I4BfBPADAH9CJfpNi2NsGqQVd4zyNIBkm3NULFjIcZ6GOj756Cex4+M7Qg0L72YbRfOhw6oapGSA8ue4tLKEmaWZmr8PxC9V9Xs24iS5So5CvAcmvIAo48ckzUeUF5CU5kMlEQ6kY5R1KFxzHU9y3ckUnuSTAL4C4CtUok7jI2tCqFgcoiKmt7M3towwpOFpCDlPHHtC6e+J5GaateRZw9TcFA7MHsDxxePYNLBJek1d7+8AxRyUo5gYnijT48s8DYUDoRv7N1b/PiD3NMQ9j+95XPo7oN546u7oxpbBLbEbOO3csjPyvqBxAOaMnyg5gTI852F+ZuRnYssIW4/5njyGuoey4WkIhUEl2kEl6q78/Doq0X+gEuW815gEER0goqeJ6Eki2lV5bR0RfZuInqt8z5uWmwSqccckNB9RMuI0vglClAJUrYhhZrS3tRsJzzUrVA441rVxDfBMg3IUne2dUpqPICUDyC3UsES4uOfFMy9i4dyC9Pe61O3+e0WC3j8uXagaP9bDU4aeJUyGqeIBFegkwr8CYJVK9BIAdwAoAviClVFdwC8x81XMLEyNDwH4LjNfAuC7lf83DUx5AZGTJ6H1FGU5AeWFoFURU4lTX6zhKZWyY1kb1yCaj1CabMm5iyClIYvrh+U0hIywZ5EqDcXNWShOkzQfYTBB85FKnjFCBpDeWQ0dpbHGk3wewFsA/DlP8vsAbLYzrEDcBODOys93AnhzyvJDkYrFoTJ5cslq3FUtNCBGRcysmYqYZoOKpyE7FAfIaT7CPAf/9WGJcEHz4b0nrHrKO66gZ5Gtg2KuiOn5aSyvLkvv8Y615lR8Crm5pDQfUeveFM1HmAwAKAwXjLWVDoOO0lihEr0dwDsBfKPyms1cBgP4eyJ6jIhuqbw2ysxHAaDyfaNF+dqI2tBFRUzSQ0uRkydhVzLVxQaoV8QIpXF6+TRePPNirHE1M1RIGwOrjjSrofw0H2FKRtwjDU/FPKgaFJ5iMA7OHQwcB1C7ASf2mHWMn4QeTZgMQfNh4lnCICtqsAEdpfEuAK8GcCtP8n4qURHA5+0MCwBwLTO/AsAbALyHiF6reiMR3UJEu4ho18yM3TfQiyiLI6wiRkdGFHRpPoJkhC2EbcPblCti/HFqHaV51zN34dfv/nXl67MKFf4tWSIckNN8hIWOGFyTNwsLZwk53nFFJcI39G1Af2d/4OcoM55UN2e/p5GE5kPL+EkhnGs9PJVS2a2y0uBJ/jFP8n/gSf5i5f/7eZL/xNbAmPlI5fsJAPcAuAbAcSLaDACV7ycC7r2dmXcy886RkRFbQ6yXqxp3TCG2CSQvvwtbbDoVMdWOdDEm9Q+mf4Cv7vkqFpcXle/JIlT4t/zKoFoREyexfWq/0vVAfVFDVE4jKm8W5GkA0XOyhn8rIc2HivEzNjSWmKkhylgEDITaVGSkVHarrDSoRNdSib5NJfoplWgflWg/lWifjUERUT8RDYqfAfwKgGcA3IdyeAyV7/fakG8TqUyehBaHiuUEaFbEIF5FjFj4zU52KJ5DRvNRvcYXRgqqiFEqofWFm6KUhreoIap6Styjw8K7dXBr5ElyQJ1KxRQ62jowMTyRiOZD1QuYPTuL2bOz1mSY7EkeBp3w1B0AbgPwCwB+HsDOyncbGAXwfSL6EYBHAHyTmR8A8CcAriei5wBcX/l/ZqC6oSeh+UjD01Cx0IB0KmJapYGTeA4ZzYf3Gv/8EfkpL8JyFJsHN6OrvauuGkonRBOVCAfkJ8mr45Osg/a2dkwMTyh5GnVd/+LOY4XwFCB/j3XlqHoBNo2fga4BjPSNWF8rKof7BOZ4kr9lbSQeMPM+AD8nef1FAK9PYwxxoGsN/Nymukc0At3GN37oLDZREdPV3hV4XZKKmFZp4FRzKn52P7blttVdE9TG9dv7vl3eiCtzK8xzaKO2uqSrSiIcKL/Hr9322shEOFDeBEVRw4a+DXXPqtOT3AuvwhHvUex5rGH83L/3/lgyVOF9j6/adJX2/SoGKZBO2a2Op/GPFa6pV1OJXiG+rI2sCZFG3FF58iQIg+ksNpWKGG84RbcQoFUaOKm0MQ3KBchoPnROOUeFpwTNR9XTiEiEA+GhkDBCRZ1EeE9Hj3LeTAZV4ycpzYdXoYfJAJJ5TVEygHQo0nWUxitRDkn9MYA/q3z9qY1BNSvSqHBQnjwGLA5lBRgVp/Z1pNOpiGml8NT63vVSmg/vNSpVRyqJ7brwVMic6e7oxtahrdVxRSXCa8YlozkJmKOFXAEzSzOh9Ph+hZOk2lDH+AHih45UwlNJaT5UZADlZwnLm5mAcniKJ/mXrI3iIkJSmg8dT+OePfdgdW0V7W16TPbKiXDFmHNN7+t8uSLm6MJRbB3aGj0WcZK6ycNTzFzenAe3Br5fsjCS10J95dgrASiU0OYu0HwMdg9GKhlxj3iPlRLhEf0+who4Tc1O4YqNV0j/rqyB0/cPfj907EnhfY/jcEOpeBpJaT5UZAAXmBoOLxzGxPBELFlR0MlpgEr0RgBXAOgRr/Ek/5HpQTUrVDb0pDQfOm6qqIgZHx7Xk6FooY0NjaGjrUOvIsajaFSUhtfTUF04WUT1gGNIyXVQIhyotYIjq6E8m+CVo1dGhrPEPf+w/x8AqCXCq0UNAYSKYa1i98/uD1QasgZOX3zmi1hZXUFnu95ZYuXwlIHzDar5hp+++NPYMlTgfRZbSkOn5PYvAfwGgN8HQADeBqA+m3cRI0txxyQxVNXFFqsiRrOMUiiN+XPzOHX2lNI9WcQaKmXHIdamTBlUK2K84SmNxLbK9eKew/OHsby6rJQIF/eY6Enuv9ffwCkuzYeq8ZOU5kPHM49L86EanjLVWCoMOjmN1/Ak/xaAUzzJJZRPh+uZsC0O3SR1rMmjKENmoepCWQFqHNjSndTeBdnMZzW8BxyPLBzB2fNnA6/xQ5bYVimh9W5QKkpDFDUo5wICvKYgpbGxfyP6OvusNXCq+1uKxg8RJc6dqK6VpZUlnFiUnkk2IkOnd0lc6CgNUVqwRCXaAmAFZaZbBw9U445LK0uxaD5UvRlB8xFrsWkoMxVyRK+lKxrfqE7qVmng5KdSmZqdCrzGD79ijspRrO9dj4GuAeyf3a+UnwBqq6Giwl/ecU3NTdUVNQRZxSqbcxBpYxLL2bb3n4YXoCrDX9RgAzpK4xuV/hkfA/A4gAMA7rIwpqaF6QSyVIaip6FD81EnQ9FCA8qL+sTiiVCaD9lGoNP1L9+TB9DcFVT+5kIypclgtEmWpIzmI2xT9yZdVfITQO3mrJIDEeMSRQ1e6PYk98K/OSah+dA1ftJIUgPxjR/VfJ7t8LcO99T/4Eme5Un+Csq5jJfyJP+htZE1IXQO4ADxXW7lyRMz4a4angDUyhVlfSJ0+ovne/PI9eSa2tOo0sOHWM5Bm20hV6ij+VBJbAuvAYg2ALw0Hyo5ECFD9iyh3FgRoVn/e5CE5kPX+ElC82E7ZKylAC0f8IucGVSiX658f6v4AvBGAK+v/OxQgeqGnkayCohvceguNkCf8lu18U2VtyqlrmS6+J8P/U/cdNdNkdeJ59gyuKWO5qPmGll4SkLzoVRC6wlPRV3vLWrQCU8BchbesK5/YUUNgV3/EngaOmsy1oauGGFIQvOhGp4CLhQ1eOnxTULF0/g/Kt9/VfL1JiujalKoehqyihjTMgAoN74JgpYCjKqI8VF+q1bEVBPIKXUl08UzJ57Bd/Z9J9ISFM8ho/nwX+OHf3OOOqcBlD+X08unMbNYPkmu6jlUw1MqeTNB8yEhVIzs+hfWwElywNF6mDVJwl2jFDyJ968TnlJhaoiLyHMaPMmTle/vsjKCFoPqBxu3WkM3PCUmz0vWvURdhoYrPNo/it6O3nBPw2cZexXN9vz20L8vNqDCcAH3P3d/5s5qrPFatSJmdGA09Dpv2XFQeEq2wflpPlQ8AbEJ7jtVJqJWDTV+7Sdfw2snXqukZERRQ114KsQq9npNV2+5uu73Qa1iBc1Hb2dv5Lj80Mo3pOAF7Dqyy64Mz7Ncsv4SbVlRiFQaVKL3h/2eJ/k2c8Npbqi6qUD5g3386OP6MjQ9DaC8OWspDQ0LLVZFjEYhgLfq6Oz5szh2+hg2D6bdZTgY4r06MHsgUml4+0Q8duQx6TWyzdpfEaN0wruycTx/6nkAip5GroiZpRksLC9obVBJe5J7IdscvaEjnRPbOsZPEpoPXS/gq3u+GoupQUcGYK/aUCU8NVj52gng3QC2Vr5+D8DlVkbVxNDZ0L0VMVoyNDwNQN960okFCzk6B7bGh8eVK2K8/cWB7FVQqXJj+U/FC5qPoGv88JY2KyXCK+/X8yfLSkMnP7Xv1D4lJSPk6ISncj25clGDhQZOfugaP0nyZjpKVtB86EBHAW4Z3FIuarC0ViJnBk9yqXKYbwOAV/Akf4An+QMArgYwZmVUTQjtjdbX+EZZjoY3o9r4JkiGjgIMq4jxW8YdbR0YHx6P18ApYxVUqtxYfv4tIKD3d0gCWeeE92D3INb3rse+2XJ4Siex/fyp57WUxvT8dE1RQ5TlHXa2Jyg8Beh/9rGMnwS5EyUZcZ9FIzylytQQFzrnNCYAeDOqywAKRkfTxNDeaBPEUNOaPDoKMKoiJm5yUygcEyfcVfD1Z7+O3/367ypfr+xpeBRnUIgmqlRVFDUoVzfli1VPQ6eEdmp2SmujXeO1mqSrCjdWIGmjROHEpfmIY/zEofnQTYQD8bx/nVxe0rbSYdBRGn8D4BEq0UeoRJMAfgjgc1ZG1cSwGXfUtZyA8galu9HqLpooL0DK3qpYDiw2oL7OPoz2j2ovtj0ze7Sqx76z7zu4/fHbpTQfQeMD1Jh+ozrSRVUdiaIG1eqmYq6olQgf6RtBX2cfVnlVy9MAap9FhRsraHOWvQdJaT501mQcmg8dL6Ba1GDR0wCSHVaMgs7hvlsB/DaAUwBmAbyLJ/mPrYyqCaG70cbhiNG1nIB4kyeOhQYEewFBcerji8cjG9/UJJA1wweLy4u46q+uwh2P36F8j3h2Gc1H2PVK4anKc2zo24D+zn75+YaoqqNT+5US4cCF3Amg5mmIuD4Qw2PWaS+bu1DU4EeSBk5+xDZ+4nj/ioqpq70LY0NjVmUA5ffrhaUXQnuXxIWOpwGe5McAfBHAPQBepBLZ4d5tQuhutEk4YnTd1CiaDz/ixIKB4MUm20RUG9/4E8g6G8fZ82exvLqMH8/8WPke3aZP4vqoxjfe5xDVYP5TzlGJcDEuZZqPyucCqCkN7z2q1wuaj7jU7X4EWdRpGj+2lVMcr0lXxjuveieefvfT6O3QL1GOgg41+q9RiZ4DsB/Ag5XvqfQMrxkH0Y1E9CwR7SWiD6UtPwhxQke6CyEOK26SrmSqi61aERNyYCtuRUxNArnSlez82nmlccXp+ifuUX2/xGcSVRGjQqUSFnaqVsRoEgoKqH6WheECAHWlIWg+/ISKUYlwQLOBUwyaD901GZepQTt0lI8m+UwqY8vgFrxs48u0y3pVoONp/A8ArwLwU57kIoDrAPyz8RGFgIjaAfx/AN6Acrnv24koE2W/sUJHmsmquDIAvYWgUw1SlaPZJ0LVqvMfilvlVUzPTyuNKY7S0O0UqMrCK6NS8VechSkDb1GDLjcUoO9paCddNXqSh23OQQonjvGju176u/ox0jcSKweoayzq0nxk6VCrjtJY4Ul+EUAblaiNJ/kfAVxlZ1iBuAbAXmbex8zLKLPs3mRD0PseeB9ue9juuUVdmo+43gxgP+Eelm+QWdCqFTH+/uKA+rN4D96pemlxwlMb+jZE3iNj+j29fLqacwCiD+1p03wMb6v+rFugoapkxD06TaJ6O3uxaWCTtqcB2C+5jlN2GydJrUvzoSvDJnSUxiyVaADAQwD+lkr0/wJQixOYw1YAXsKi6cprNSCiW4hoFxHtmpmZiSXonw/9Mx7Y+4Dy9XErm3QmTxxPQ6XxjQk5heFC4OYs20SIKJCDyYskDZyEAlhaWcLMkto8qCa2NZpEid4lkUy/vt7XgH4CWScRLujxgRiehuYmeHzxOJZWlgCoHz6M1fUvhteoHTKOoZh0DSxAPwzWjJ7GTQCWALwPwAMAnkeZtDBNyN61ul2KmW9n5p3MvHNkZCSWIF2LI25lE6BvPelMnjjlinE9jbCKmKCe0Tr9xSeGJ9BGbcZDR7J7dK7v7eyNLGrwK07ZxqFSqjqzNBOrJNbW9UA9Q6wSN1bA+gqyqPO9eQx3D9sP51byZjpMDdpVWgm8/yxAaWZQidoB3MuTvMaTfJ4n+U6e5I9XwlVpYhq1LWbHAOgdqVaELs1H3I0W0LBqY06c2NZTHAWoyd6qc76hs71Tq1zR+37peicymo+g66vU7RFUKlFtTFWrjgD9kljV64d7hpHvyeuFp3zlwEC00ikMF3Bo7lBNUUPUvWkZP7o0H3GS1Lo0H00XnuJJXkW5zeuw5fFE4VEAlxBRkYi6ANwM4D4bgnRpPuJYNbo0H3FkANE0H0FytGSExJzDQg5RFTHSqiNNBRA0Lhm8z65zYj3KM/U/h6D58JIQAuEbnLcayqbnUMwX4+XNPJ0CVfqLr/IqDs1diDZHdRnMqvevm6Rub2vHtlx0aDaJDJvQCU+dBfA0legOKtHHxZetgcnAzOcBvBfA3wHYA+BuZt5tQ5aJ3sRR0KX5iGM5AYhsfGNCTli+IbK5kE7VkUa5oq4CEGOt3qN4Yl0ciguriAkqO/Y2VgLCN/caT0Mzsa3zWb5k3UvQ3d6tfH21qEHD05B5plH9zOPSfOggbrVhLENOM9TWVJ5GBd8E8IcoJ8If83ylCma+n5kvZeYdzHyrLTnalToJNnSdpCug72moNEpKKqdK8yGRERSrV0luyhLIRxaOKNF81CgADaUx2DWofI+XhTesqCGQf8vTWAkIf88FzQdg77AeAPzxL/8xPvvmzypf782bqbaXFePyHwoMG6suzUecNRmX5kN73We0E6UKIvtpCPAk32lzIFmDv/FNFJKEju599l7t8enKAMoLVNb4xg/TCnCN19BO9YeMZBtH3Vj8CeTKs0zNTuGyDZeFjse7geko/5H+ETBYixsrqvFNEP/Wvc/eizVeU2rLKjya3TO7lZXAlaNXopgr4rL14e+VFzvW7cCOdTuUrwcubIKq7WXHh8brihqq8y6KSmV2f2jvkurfi7Em49B8xDp0m79A8zHQNaAko2nCU1Sir1OJfpVK1Cn53XYq0R9RiX7bzvAaB12aj9gbbU6d5iPJZg7EKPGLoQBlCiAoEa7S+EZ2vgFQzzcAZaoL1YqYqhJQtARFriLKm5PldQq5ApZXl3F04ah6ArkiR/WzWde7Dvv+YB9+fuvPK10fF9VyYIUwG1AuahgfGpeGp+K2ivUjDe8fiB+eAvS8/2YKT/07AL8I4CdUokepRPdTif6BSrQfwF8BeIwn+TNWR9kgFHIFfTdVd/IoWNsCcb2ZKJqPIDm6CKL5CEqEqzS+kSXCAfV8A1C2nFUrYqrhJsWNQyjEalGDDv+WRwFGxfOr98RIbKeBYr6IuXNzOHnmJAD1Fqs6SiMuzUeccK7NRLiQAejR1TSNp8GTfIwn+T/zJO8A8DaU6UTeD+AKnuTreZLtxlYaCJ24Y5KNFtCznONMHq1nSWChyWg+IpsLhZ1v8Hkpmwc3o7u9W9kLAIDtue0A1BWN8DRUkq7i2aKKGqSJcI8CVLXQ4+Qo0oC/J7kqN5b/JDkQPO/6u/qxsX9jrEIIHYi8mQ7NR1xjUWd/aSZPwwsG0M+T/CQAphINmh9SdiAqYlRoPhKHjhRbnwL6E1TISSM/A0iaC0Wwt4Ztzn6F00ZtyuWK4jlEfF5ls/Geu/DTfATJ8HJjhYWn/HNjW24bAGjlAuJUQ6WBak9ynaZPuSKOnj5apcdXeQ9SMX40aT7iKCdR1JC1TpQq0GG5/XcAvoxySAooH6z7moUxZQbFvPrkibvRVieP5UqKOOWKcV1uaRvToORmREVM0gZOYlyqRQ3ecxeAJqFiSBtT2XP0dPRgy+CWWFVHWfc0dKjbp+amAEQf7hP3WD/XFOPQre5a0e1J3lThKQ/eA+BaAPMAwJP8HICNNgaVFejEz+NaNTo0H0nDU0E0H0FydBFE8xF20jkqTh2UQNZ5v7rbu5UrYrznLsLGVXN9ZVMSNB+yxjdRzYVUE+E78jsw2DWI0f7o6qE0IWg+nj+l52kA0Co71qH5SOJpeMcVKSdm6EjX+2/G8NQ5nuRqnIZK1AEJ71MrIU7VUazJo2g5Jw1PAepnD+LICaL5CO19HWHRBzVwOnnmJObPzYeOxxvuULVQxeaueralhoU3pKghlH9LIxE+2D2Iqf84hZtfdnPUo6SOQq5QVRoqG7XfYFAJTxVyBW2aD915rEvzEdcL8BsMNmTYgI7SeJBK9F8B9FKJrgfwvwF83c6wsgEdmo+4STdAneYjqacB2PWahBxpeCrgb0V5GqENnBRCR4BeKECEkfw0H2Eyqu1oQ97jMP6t6flpnFstJ11VLPR8b95Kc52kKOaLWonwalFD5f1SoSDRmscx12Qcmo+4xuLC8kK14syGDBvQURofAjAD4GkAvwvgfgD/j41BZQU6NB+JNlpNmo84k0enXDFxwl0jET7QNYCRvhGt8w06rWIBVD0HlcY3NeEmhfCBv0kUEEylEqQ01nit2pM8a7kKHYgzR4Dac/iLGpQS4THO6STxAlQQu0orhvefBSjP0ArD7ad4kt8G4BYAP+TJDPH1WoJyvX6SjVbRekoycXo7ewNpPoIQd7F5K2IANfZWnQZOOq1iAVSroVSKGvxNn3So28MqYqL4t3QSyFlFnPayXg9QJa+jQ/ORdE2mEZ4C1L3/pgtPUYn+iUo0RCVaB+BJAH9NJbLb2i4DUM43JPQ0AAXLOYEMIefAXLgMr5xYMgTNR6UiBlBsLhQUnpJ4Ket712Oga0A5POWn+Yi6x6s0puamajisZNeLZwsLg0Xxb+mEdbKKWO1lPRVnKolwHZqPpGtS0HxEykmQCAfUPY2sGBQ6M3SYJ3kewFsB/DVP8tUo9wlvaYRVxMiQyNOIan2awHISctJKuHsVYGQb05CKGJnCUc1ReC1XVauuJkeRL2J5dTmUHl/GwhvI9Ct5P8eGxtDR1qFVdZRVxKJuz18oalA+q6JJ82HT+wfiewFD3UNY17uudT0NAB1Uos0A/i8A37A0nsxB2QtIEDoSjW9UKnWABJ5GAM2HSTmBbUzDkpv5Is6vnZdWxAQmkBXzDUD5OVQrYmoO6ylsHLLmSrKKmKDnEHmzangqIxtDHIi8GaBP3a51Kl7T+IkDXcqSJIacivefRIZp6CiNP0K5j8VenuRHqUTbATxnZ1jZgU7pJZAsdGR9ggbQfPiRxNOQ0XyotDEF9Bs4RZUrejch1YoYfyIc0KNuL+QK0oqYqLMqreBpCJoPQJ+6XZd/S4Xmw0TIOA3lZFuGaegkwv83T/KVPMn/vvL/fTzJv25vaNmAbqVOotCRYngqLnQPLcVZbDKaD9U2ploNnHJFLK4s4oWlFwL/rj/cEXZiWzbWbcPbQKBY1O06ZceCrgbIjjUZF1WaE41EOKDXKlaVqSHJmhRFDUpEoglCR2JOhuXNksowDZ1EeA+V6D1Uok9QiT4jvmwOLgvY2L9RieYjqadRyBUiaT5MeDOAQu4kYVGcP3wQlQgPq4gJTCAregHAhU1DqRrKE0bq7uiu0nyEyahh4Q0IZ0bxb6mGZrIOXZqTdb3rMNg1qEelopoDTLBedM/2JPH+z62ei2RqaNZE+N8A2ATgBgAPosw9tWBjUFmCDs0HkMzTiKL5SOrNyBrf2JDjf7+iEuFd7V2BvUsCT1IreE1+y7WYjy5qkIWbtNrRBowrlH8rRtVRVqFL3U5EdafiVZl+lT3mBBu6apVWEk8DUAt/N52nAeAlPMl/CGCx0sXvjQB+1vSAiOgjRHSYiJ6sfP0rz+8+TER7iehZIrrBtOwgqNbrJ5KhYDkn9TRkjW9syPHTfEQlwsU90lLVAIWjkqj0x8hVQo114aaIjcMfeqsWNehQqXjPN2RkY4iLOCy8QjGrelubBzYrFzUkgQ7NRxLFBKgl3JvR01ipfJ+lEr0MwDCAgvERlfG/mPmqytf9AEBElwO4GcAVAG4E8AkiSQ9RC1Ch+TCx0QJq1lOSyaNiPSXOz/gswahEeHVcGolwQfMRpgDqchoKFqqsGmp6fhorqyuB18uaK8lyGlH8W96xNiuEMtd5DrG+RMl11BpSLWowsSZVaD6MVGlZNkpNQmeG3k4lygP4QwD3AfgxgI9aGZUcNwG4i5nPMfN+AHsBXJOGYBWaDxMhHSDC0zAwcXToERIrQA89RFvEVAuqiIkK66i8X2IDU/VO/J0C13gtMOmqSt0e9hyj/aPo7eitGWuz4srRKzHUPYRtw9uU7/HT48dp4CSDMeNHQTnFXSs9HT3YPLDZqgzT6FC9kCf505UfHwSw3c5wqngvEf0WgF0APsDMpwBsBfADzzXTldesw+sFrOtdJ70mqVXT29mLTQOblGL0SSZPIVeo0nz0dvaGyokLv0WvGp5iMKbmpnDp+ksvjCUigfzksScD/6Y/sarS+EZ2WA8obxyimZNfRh3NSa6Ib/z0GzXeRdhziLzZnhf2ZCYEERejA6OY/S+zWnO02sBJk1b9qz/5aug1Jr3/nVt2BstJmKRW9f6zMjeUlQaVqBvAr6Mckqrex5P8R7pCieg7KCfV/fhvAD6JcktZrnz/MwC/DUjfMenuRkS3oMyPhYmJCd3h1cG7cVy95erEfy9QTsRBn6SWk5ABlGk+XrrhpVbkVGk+PJxCOslNr9IIDevkirj32XsDr/EnwpV7klOtpyHGJUMQC6+oiNkyuEX6d+ueJV/Enhf2NL2nAcTvX6HbwEnQfAx0DegPUmVcOvmGBIZcMVfE9w9+P/SaLHkaOjP0XpRDROcBLHq+tMHM1zHzyyRf9zLzcWZeZeY1AJ/ChRDUNIBxz58ZAyDld2Dm25l5JzPvHBkZiTPEGihV6pjY0CN6PpjwNFROuCeVIzZnL6eQahtT/7ii+ouH0XzIqnF0E9uC5kOnSVTQqXiVA45Z2RjSRFxPA7C7JlVpPkyUqB+aPxSYNxNoOk8DwBhP8o3WRlIBEW1m5qOV/74FwDOVn+8D8AUiug3AFgCXAHjE9niA4IoYL4xs6LkivvTMl3B+7Tw62oI/GhOeRhoJd2E5Rp3TAIIb36hUHe0/tR9jQ2N1v5edMC7minjwwIOBlps/jBRFjx+UCAfKFuq1E9dGPof3WVrB09DFQNcANvRt0KJS8b7HPzsqL+I0Fc61HToq5st5s0Pzh7A9L4/8N2si/F+oRMZLbCX4KBE9TURPAfglAO8DAGbeDeBulBPwDwB4DzNH93w0BNWka9INPYzmw8TEkdF82JDjLVdUqZ4KqohROd8QxpALoK4aKqwiJqhTYGB4SvJssoqYSKbfjPb+TgvFnF4Dp9S8f0VizKTGIhDt/WfFC430NKhET6OcO+gA8C4q0T4A51DOMTBP8pUmB8TM7wj53a0AbjUpTxWFXAG7T+wO/L3J0NH+U/tryN9MypDRfNiQ46X5UEmEi3t0GjgJmo+wBk5A7abhDc+t71tfd480R5Er4r6f3hcow3+9qIipYfqNUJxXjFwBAiHfkw+8ppVRzBfx6JFHAagpjQ19G9Df2W81zArIixrq5CT0NGqMjKL8miwlwlXMmjcB+FUAbwDwEgC/Uvm/eP2igOitkDR+GYbI1qcGLCcgulzRiKfh8QJUEuHVcWmcb4ii+ZDmNCIoKGTyCrkCTiyewNLKkvR62abkD2tE5TQu23AZnvv95/C6wusCr2ll6DZw8p4ktwlVmo8kiml8eBzt1B5pyGXF01BRGsdRzi38J5QP1R3mSZ4SX1ZHlyFE0XyY2NCjaD5MWE5AtMvNnNyq8VpPKolwcY+/8Y1K1VFYf3GgPhEuxiVDWKdAmVWrSt2u4m3tWLcjMxtD2ojViyNqHhusNlQpUImLjrYOjA9HMzVkBSqfzp0AdqLcG/wNKJfAXnSIjJ8b2NBVaT6Swtv4JggmFBNwgfJatYwSqF2gKgnkyPCU51mqFTEanQKjqNuDuLEOzR2q9i5R9bYuVnjDsVpKI4Tmw2jIOEI5peH9N1N46nKe5N/kSf4rAP8ngF+0PKZMIrJe31ToKMxyNhieAqKfJQkEzYfgFNJKbvot9JDnLeQKOLxwGMury3W/C+IyCrNQo6qhZDKC+n2s8ioOzR1Seo6LHV4qFdVNPqh3iYCJ9aJC82EidKTk/WfEC1VRGtXiYZ7k8JZvLQzVTl5GJo/t8JSC12Rig/Oyl2qVUXo9jQiFI2g+xObsRRBrahhzrSyUJmg+tJpE+d5jVcV5sUIUNQDxGjiFIcl6UaH5MOJp5Is4dvoYzqycsSbDFFQ+nZ+jEs1XvhYAXCl+phIFxzdaDFE0H6YS5MVcsUrzUSfDkKcRZT0xzFg1XqJHlY2gSvNhsIETUP9+hTW+kY01jB4/rEkUUEul4pRGMERRAxCvgZMMxtakQsLdVDg3qBosS55GZMktT3IqTLLNgDAX0mR4CgAOzh3EZRsuq5VhyNPw03z4YczTyBXxtZ98DYCa9Sij+YhMhIdsHEGd4GQ0H1HygjaOIIXor4hR9bYuZhTzRRxeOGzM0zAZzg2j+TChnLzP8jMjP5P479mEM300oFKpYzKBbAt+mo+ga5KimC9iZW0FK2sr6tZjvnZcUV5KGM1HUM/pUEUTkqPQSYT7K2JcIjwauqfio2g+TK7JMJoPU4lwIML7b6LwlEMF/ooYL0x7GkEbmgkZQk6UhZZYRtwySk9FTFQCOYzmIyinERXSClIac+fmcOrMqZrXo1h4a5h+M7Los4o4/FtpeP+FXKFK8yGVYyB0tGlgE3o6esK9/4x4qk5paMBfESODickTRPNhynICwssVTSbCBXSUhrciRpXsMEzJqtB8CASW0AYomigWXpcIV0ccKpVU8g0KzbuSrheRNwvMaThPozkRWnppKOnWRm2BSVejnoaH5kMmx4Ri8lbEqP69OFVHQdZmUCJcRvMhEHhYTxI+iFLioiJGnCR3SiMc12y9BgNdA1LyySAUc0VMzU4FFjWYQFTI2KRnHqYAnafRhBAWqnSzMR06CrNqDOUbgOCFYOI5vBUxOp4GgBovKLKBU76IE4snsLhcy9QflAgX9+hUQ8lKrqN6WvsrYrKy6LOKy0cux8KHFwKZXmUo5oJpPkytSVHUYLuyKawU3CZ9kS6c0tBAGM2HydBRYTjA0zA4cVSqjozIqSgn1YXr3ZyjNuWqjJBeHEH3h3knsta0+d48hruHa2SE/X3gwrM/f1K9T4SDHkJzgIbWZBTNh6nQUTFXxKmzpzB3ds6aDBNws1gDYTQfpj0NGc2HaRlA8ClnU1axUAKqG2a1d0mFs0rl3iBuqKDqKSC4qCE0sS3hkwKCPw/x7DrNhRz0EBY6Mh3ODfMCbHv/LhHexIhKvJlKUgP11pNJb0Y0vgmy0ExZNXGaC3lPkgPRiz5o4wjzVIp5eVFDaKdA38YRFv4CLhQ16LQxddDDttw2ABFJaoMHVQNlGFJMQHDlZFbmj1Mamgis1LEQ0gnaBE1NnrCEuymrJkkZZdSmLLCxf6OU5iNM6QTRwkR1Cjwwe6CmHBgIfjZR1OA8DXsIo/kwvSaDaD6MJcKj8ozO02hOBNF8mHaFgWDryeSGbjMRDsQso6xsziJ0FPW8QTQfUTkNoP49Di2hzRdx5vwZHF88DiA6ES7u0Wlj6qCPwNP6FtZkUMWdic8235PHUPeQdaM0KZzS0ITYBKfmpmpeNxk6Wte7DoNdg1Ir2CREueLqWm3XXBsJd52FW8wXsby6jCMLRwAotv+UbBxhnkpQ45vQ9rISPqmo8Xk9U+dp2EGU928932AodCSj0TEtwwTcLNZEYL7BoFUjupL5rRrT4SlB8yE2Z68cU1bx+PA4/u3L/y1ev/316uOqvMe6PaMDw1OSZwmqiIlKhAMXNg6VnIsoCVV9Dgd9BNF8pOX9m0xSh3GcZcVTbcgsJqK3EdFuIlojop2+332YiPYS0bNEdIPn9auJ6OnK7z5ODXoH04o7yiwOk5aTkAHUu9wmE+Ft1IZP/dqnsHPLzuiLxbjytUpDqYGThOZD5RyFTNEEvb/+k+QqOZeaPhEZsRRbDcV8MZTmw8R6GR0YDWRqAAwacr68mYDzNIBnALwVwEPeF4nocgA3A7gC5daynyAiwbL7SQC3ALik8nVjaqP1oErzEVDZZAoymg8bngYgT7g30qqJU6oqe5bIcxSaPcn7Ovsw2j9a72kosPCGjcMhGaKqDU1AhanBBIq5IpZWljCzNFP3u4va02DmPcz8rORXNwG4i5nPMfN+AHsBXENEmwEMMfPDXJ4JnwPw5vRGfAFBk8fGhh5E82Fq8giaD5nL3UirRlTE6CSQZRtHVPhIVhETxUbr/exVE+ECWVn0rYa0qg2DmBpMh6cAuwowKbJm+mwF4PUxpyuvba387H9dCiK6hYh2EdGumZl6jZ0U0nyDpdBRDW2F4YkjaD5sJ9zjoJgvJvY0osJH4j32FjVEsdF6Nw6VRLioiFF9Dgd9jA2NSWk+bKzJIAoh0+eapN5/q4eniOg7RPSM5OumsNskr3HI61Iw8+3MvJOZd46MjOgOPRIymg8bZyiA+tanJmUAAVVHDQ5PAeXn16HfyPXkkOvJyT2NEEJBQL+97MG5gzi/dl4pES4qYlSfw0EfHW0dUnp8455GAM2HSU8j6LDiRZEIZ+brmPllkq97Q26bBjDu+f8YgCOV18ckrzcEQTQfpmUAcsvZ5OSRJYNNJsLjQrDwAnrtPw/MHaj+XyWnAdTnQaKUxiqv4vD8YeXDh7r8Ww76iCL5NCUDsNsgbaBrACN9IxenpxET9wG4mYi6iaiIcsL7EWY+CmCBiF5VqZr6LQBhyscq4tBk60JG82HF08gVMT0/jeXV5Ro5jbZqYjVw8m0cUTmHauMbn3cSmtj2bBwqiXCg7JmGjcMhOdKsNpSV29v2/oHs5MQaVXL7FiKaBvBqAN8kor8DAGbeDeBuAD8G8ACA9zCzOHn2bgCfRjk5/jyAb6U+8AqkXoClDd325CnkCmAwDs4dNPY3TSBOA6fCcEFO8xHwmRARtg1vq/PmojwNALEIFbOy6FsRhVyhvqghrWpDw6GjIO8/K2hU9dQ9zDzGzN3MPMrMN3h+dysz72Dmy5j5W57Xd1XCWzuY+b3cwHcxDU8DqLc4bDyyNK6fkfCUgE4DpxqaD8VzFP7wVNizTwxPlOnxY1C3O0/DHmRnjkyvySCaD+OeRiVv5mVqcOGpJoeM5sOWp+Gl+bAlA6j3mhptFQuaDyBeAydAn+ZD3BN2fWd7J8aGxrRYeEVTofbqkSMH00jD+w+j+TBtLK6sreDwwuHqaxdFIryVIWg+bCbEgPKGtrK2gqOnjwKw482MDY2ho60jc56GoPkANBLhQTQfEYfvvBUxyu1lPQcvo66/fORyfOz6j+FNl75J6Tkc9BFF8mlMTgjHmTEZskiG8zSaH0G9FUxbHICHtsKCp9He1i4tV8yCVaMb1qmj+VB4v7yKRqe9rNfTiBofEeGDr/kgRvrNl387lFEtarBcbejPmwF2EuGA3SqtJHBKIyb8PR/SCh3ZgP/QkklahCTQVRpJaD72n1LPURSGCziycARnzp+J/PsO6UBGj29jHhfzZZqPE4snLsgx7JlPDE+AQHX5mazMM6c0YkJMHsERY8OqEZPHT5BnevL447RZCE8B8aqOvOGDqPyEV8aB2QPa1VCO8jxbkHn/pudxGjnArvauat6sRkYG1iTglEZsBDHEmvxguzu6sXVoax3XkfGFkC/ixOIJLC4vVuVkwaqJ1SrWs3FElc8CnoqYGO1lXRvXbEFG82HcwArghrKxJrOWZxRwSiMmgvINplFDkGfR0wAuKMCsTFCRo9Bq4OSj+Yi611sRo3vC27VxzRaK+fqiBtMQc1JW2msSdd5/Rgw5wCmN2PD3mLYaOrKYCAfkibcsTNBXbH4F/vXP/mtcO3Gt8j3FvIfmQ6ESStzjPawX9exbBregq73LKY2MwR86smH8yGg+bGzohVwBh+cP49z5c0b/rgm42R4TfpoPaxu6hObDNPwlfllJhPd29uLzb/18VUGrwLtxqOQ0xD061VBt1IZtw9tc7++MQeb92/hsZIdubax7L1NDVrx/wCmNRPC6kNY8jfyFyWNLxsb+jejt6K3JnWRlgurCu3FE8UgJFHIFLK0sVU+S63gnqtc72IfM+7cxj+sS7haUk9/7d+GpFoHsoI+1ag1PSahpGf5yxSyV9+lifGj8As2HQiIciJfYLuaKF0pum1TBthpkNB9WPA0JzYfNdQ84T6Nl4KX5sBXSiXP4LK6cGrc+IxNUF36aDyUFIBLbGv07XBvX7MFP82FzTXppPmwkwrcMbkFnW2fm8oyAUxqJIGg+jiwcsbahbx3cis62zppqDWsut2XFlBZE+EA1p5GkJ7nq9Q7pwOv92wxPAXZzJ+1t7diW22ZdAcaBm+0JUOMFWAodeWk+bE6cYq6I+XPzOHX2VKZomONAbByq1VOiIiZOT3LV6x3SgTirwczW8gCNyJ248FSLoCbfkELoyLYMoFx/3szhKaD8uRxZOIIzK2e0aNV1PA1vRZfzNLKDYu4CzYetjdbP1ADYy51Uz065RHhroDp5vHFHm6EjS96MkAFcUIBZmaBx4C271aFV10mEb+jbgP7OfgBOaWQJ/qojG/M4iKnBNIr5ImaWZnB6+bTzNFoF3sljO3Tkpfmw6WlUF1tGJmgceE9s6ygN0fNd5R5Bjw8093vVapBVG9qSk1buRHgbWTHknNJICG8MFbDzwfpjqDYmaK4nh1xPzvpiSwNisR2cO6jdiwPQb/rkPI3swEvzYdM6r6s2tGnIZWxNNqpH+NuIaDcRrRHRTs/rBSI6Q0RPVr7+0vO7q4noaSLaS0Qfp4yo3Wq+wWboKIAkzbgcTwVVRt7eWNg8uBnd7d3K1VNAzPaylXua+b1qNfR39WNj/8aq92/rsxF5s3Pnz1lTTl5j0YWngGcAvBXAQ5LfPc/MV1W+fs/z+icB3ALgksrXjfaHGY3CcKGG5sNWQgwA9s3apa3wVh1lZYLGQRu1YVtuGwC9/uLe+3XucZ5GtlBj/NjyNHw0HzbW5EjfCPo6+6pGaVaMk4bMdmbew8zPql5PRJsBDDHzw1yOA30OwJttjU8HguZjanYKgB1PY2P/RvR19lmn4hZdyVTpN7IM3dDR+NB49X3V9U6aWcG2Igq5QtUrt2lgAbCaz/QzMGdlnmXRRCoS0RNE9CAR/WLlta0Apj3XTFdek4KIbiGiXUS0a2ZmxuZYU+muJ2g+bBPkFfNFnD1/FsdOH8vMBI0LXaUhihoAdSXw+u2vx/te9T68auxV8QbpYAVeenybMgBPtaHN3InFSrA4sKY0iOg7RPSM5OumkNuOAphg5pcDeD+ALxDRECD9RALVOzPfzsw7mXnnyIjdvszC4rC+oVfqzwF7lq2soX2zIm4vDkBd0Qx0DeC2G25Df1e/9vgc7MFL82FrrXhpPmznTrKWCO+w9YeZ+boY95wDcK7y82NE9DyAS1H2LMY8l44BOGJinEkhaD5sh468iVpbEApwen4al224zLo8m4iTbyjmi/jewe9lxqJziAcvAaWtz7KGqcFy7mRheQFU+ZcFZCo8RUQjRNRe+Xk7ygnvfcx8FMACEb2qUjX1WwDubeBQqxCTJ43QkYAtGcI6b/ZEOBC/VazuPQ7Zg9f7tzmPvWW3ttf9RZ8IJ6K3ENE0gFcD+CYR/V3lV68F8BQR/QjAlwH8HjOfrPzu3QA+DWAvgOcBfCvlYQeimC9bA1ZleEtCLS2Evs4+jPaPlmVkZILGRfXgncZzOKXRGhBMDWmsyTQO9gpkxZCzFp4KAzPfA+AeyetfAfCVgHt2AXiZ5aHFQhofbBqehpBzfPF4ZiZoXKzvXY+BrgHt8BSQncXpEA9d7V0YGxrDoflDdtdKrogXll5AG7Wlsu6zAmdSGUAajKdegjyrLncKuZM0IMoVdd6rl296OW7YcQOu3nK1xZE5pIE0DAAhw2aJ+lD3ENb1rgOQHe/fKQ0DqPECLE1SQfNhG610yvnS9Zeir7NP+frB7kE88JsP4NL1l1oclUMaSGMepxU6ytp5IKc0DCCt3gqpLIQWCtHcdsNt+Ju3/E2jh+HQAKSx0aYVOoqTn7MJpzQMIPXJk4ZVk5EJmgQTwxO4YuMVjR6GQwOQxpoUNB+A3fVSGC6UZWTEkHNKwwBqJk+Tb+je3ImDQ7MijbUi8mZAa3g0qnBKwwAEzYf42RbSSFJPDE9YrQZxcEgDcRgB4iCN0FHWvH+nNAwhDYvj2olrMdo/iq2DgbRbidHZ3omxobHMTFAHhzgQNB+253GankZWDLmGnNNoRaRhDVy16Soc++Axa39fYHt+OxbO2T0Y5eBgE+1t7diW21btdmkLaYRzC7kCCIT2tnbrslTglIYhZC3umAS3/vKtTmk4ND2KuSKeOfGMdRmAXWOxp6MHd775Trx6/NXWZOjAKQ1DuGbrNehq78LG/o2NHkpivGb8NY0egoNDYrxy6ysxe3bWqoy0Qkfv+Ll3WP37OnA5DUP4hYlfwMKHF7BpYFOjh+Lg4ADgI6/7CP7ld/7FqoysJanTgPM0DKKrvavRQ3BwcKggjRzAcM8wPnb9x3DjSzLRfToVOKXh4ODgkAAffM0HGz2EVOHCUw4ODg4OynBKw8HBwcFBGU5pODg4ODgowykNBwcHBwdlOKXh4ODg4KAMpzQcHBwcHJThlIaDg4ODgzKc0nBwcHBwUAYxc6PHYBVENANgKubtGwC8YHA4aaPZxw+4Z8gCmn38QPM/QyPGv42ZR/wvtrzSSAIi2sXMOxs9jrho9vED7hmygGYfP9D8z5Cl8bvwlIODg4ODMpzScHBwcHBQhlMa4bi90QNIiGYfP+CeIQto9vEDzf8MmRm/y2k4ODg4OCjDeRoODg4ODspwSsPBwcHBQRlOaUhARDcS0bNEtJeIPtTo8eiCiD5DRCeI6JlGjyUuiGiciP6RiPYQ0W4i+oNGj0kHRNRDRI8Q0Y8q4y81ekxxQETtRPQEEX2j0WOJAyI6QERPE9GTRLSr0eOJAyLKEdGXiegnlfXw6oaOx+U0akFE7QB+CuB6ANMAHgXwdmb+cUMHpgEiei2A0wA+x8wva/R44oCINgPYzMyPE9EggMcAvLlZPgcqN43uZ+bTRNQJ4PsA/oCZf9DgoWmBiN4PYCeAIWZ+U6PHowsiOgBgJzM37cE+IroTwPeY+dNE1AWgj5lnGzUe52nU4xoAe5l5HzMvA7gLwE0NHpMWmPkhACcbPY4kYOajzPx45ecFAHsAbG3sqNTBZZyu/Lez8tVUFhoRjQF4I4BPN3osFyuIaAjAawHcAQDMvNxIhQE4pSHDVgCHPP+fRhNtVq0IIioAeDmAHzZ4KFqohHaeBHACwLeZuanGD+DPAfxnAGsNHkcSMIC/J6LHiOiWRg8mBrYDmAHw15Uw4aeJqL+RA3JKox4kea2pLMRWAhENAPgKgP/IzPONHo8OmHmVma8CMAbgGiJqmlAhEb0JwAlmfqzRY0mIa5n5FQDeAOA9ldBtM6EDwCsAfJKZXw5gEUBD86xOadRjGsC45/9jAI40aCwXNSq5gK8A+Ftm/mqjxxMXlXDCPwG4sbEj0cK1AH6tkhO4C8AvE9HnGzskfTDzkcr3EwDuQTn83EyYBjDt8VK/jLISaRic0qjHowAuIaJiJel0M4D7Gjymiw6VRPIdAPYw822NHo8uiGiEiHKVn3sBXAfgJw0dlAaY+cPMPMbMBZTXwD8w8282eFhaIKL+ShEFKiGdXwHQVBWFzHwMwCEiuqzy0usBNLQYpKORwrMIZj5PRO8F8HcA2gF8hpl3N3hYWiCiLwJ4HYANRDQNYJKZ72jsqLRxLYB3AHi6khcAgP/KzPc3bkha2Azgzko1XhuAu5m5KctWmxijAO4p2x/oAPAFZn6gsUOKhd8H8LcVI3YfgHc1cjCu5NbBwcHBQRkuPOXg4ODgoAynNBwcHBwclOGUhoODg4ODMpzScHBwcHBQhlMaDg4ODg7KcErDwcHBwUEZTmk4ODg4OCjj/wf+Fy4gugZLrgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "f,H = signal.freqz(h,1,512,2)\n",
    "angles = np.angle(H)\n",
    "fig, ax3 = plt.subplots()\n",
    "ax3.plot(f, angles*180/np.pi, 'g')\n",
    "ax3.set_ylabel('Phase(radians)', color='g')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9ec012cf",
   "metadata": {},
   "source": [
    "#### 很难将 360° 跳跃（由 angle 中反正切函数的定义导致）与 180° 跳跃（表示频率响应为零）区分开来。\n",
    "\n",
    "#### unwrap 消除了 360° 跳跃："
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 275,
   "id": "76986f5b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 3.0)"
      ]
     },
     "execution_count": 275,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZwAAAD4CAYAAADYU1DBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAqcklEQVR4nO3dd3gVZfr/8fdN6KjIUkRBBCWAyIJoRCysooCILoiKYgNFV7pfsfckq7iiC+wamgIqqIBiAaQIIooNgVClJ4qr9CYtIC33749z4i9CgCTknJOcfF7XNVdOppznHmbNZ2fmmWfM3REREQm1IpEuQERECgcFjoiIhIUCR0REwkKBIyIiYaHAERGRsCga6QIipUKFCl69evVIlyEiUqDMmzdvi7tXzM22hTZwqlevTnJycqTLEBEpUMzsf7ndVpfUREQkLBQ4IiISFgocEREJCwWOiIiEhQJHRETCImoCx8xamtlKM0s1syciXY+IiPxZVASOmcUAA4FrgbrAbWZWN7JViYhIZlEROEAjINXdf3L3/cAYoM2xNtiyZwvpnh6W4kREJHoCpwrwa6bf1wTn/YmZ3W9myWaW/L/t/yPu9Thm/jwzbEWKiBRm0RI4lsW8I94s5+6vu3ucu8fVKFeDrXu3cuWIK7nxvRtJ3ZYa+ipFRAqxaAmcNcCZmX6vCqw71gZ/KfUXVnRfQe+revPZT59Rd2BdHp76MNt/3x7KOkVECq1oCZy5QKyZ1TCz4kB7YMLxNipVrBRPNXmKlJ4pdGjQgf7f96fmqzUZOGcgB9MPhrxoEZHCJCoCx90PAj2AqcBy4H13X5rd7SufVJlhrYcxv/N86p9Wnx5TelB/cH2mpEwJVckiIoVOVAQOgLtPdvda7n6Ou/fOzXecX/l8Pu/wOeNuHceB9AO0GtWKlu+0ZOmmbGeXiIgcRdQETl4xM9rUacPSbkvp16Ifs9fOpv6Q+nSd2JXNaZsjXZ6ISIGlwDmK4jHF6XVJL1J7ptItrhtD5w+lZlJNXvn2FfYd3Bfp8kREChwFznGUL12epFZJ/ND1By6vdjmPTX+MuoPq8uGyD3E/oue1iIgchQInm86teC6Tbp/E1DunUrpYaW4eezNXjriSeevmRbo0EZECQYGTQy3OacGCzgsYct0Qlm9ezkVDL+LucXezbtcxH/sRESn0FDi5ULRIUTrHdSalZwqPXvooo5eMJjYpln/O/Cd7DuyJdHkiIvmSAucElC1Zlj7N+7C8+3JaxbYi/st4ag+ozTuL39HAoCIih1Hg5IGzy53N2HZj+erur6h8UmXu+vguGg9rzLe/fBvp0kRE8g0FTh5qclYTZt83m5E3jGTdrnVc/ubl3DL2Flb/tjrSpYmIRJwCJ48VsSLc1eAuVvZYScIVCUxKmcS5A8/lielPsHPfzkiXJyISMQqcEClTvAzxV8azqscqbq13K32+7UNsUiyvz3udQ+mHIl2eiEjYKXBCrMopVRhxwwjm/mMutcrXovPEzjR8rSHTf5oe6dJERMJKgRMmcWfE8dXdXzG23Vh2799N87eb8/fRf2fllpWRLk1EJCwUOGFkZtxc92aWdV9Gn2Z9mPnzTOoNrscDUx5g656tkS5PRCSkFDgRULJoSR677DFSH0jlvob3MXDuQGKTYvnP9/9h/6H9kS5PRCQkFDgRVKlMJQZfP5hFXRZxUZWL6DW1F/UG1WPCygkaGFREoo4CJx+oV6ken97xKZNun0RMkRjajGlDs7ebsWjDokiXJiKSZxQ4+YSZ0Sq2FYu7LCbp2iQWblhIw9cact+E+9iwe0OkyxMROWEKnHymWEwxejTqQWrPVB5s/CAjF40kNimWF79+kb0H9ka6PBGRXFPg5FPlSpWj3zX9WNptKVfXuJqnZzzNuQPPZcySMbq/IyIFkgInn4stH8u49uOY0WEGp5Y8lds+vI3L3riM2WtmR7o0EZEcUeAUEE1rNGXe/fMY3no4q7evpvHwxtzx0R38suOXSJcmIpItCpwCJKZIDJ0admJVj1U83eRpPlr+EbUH1ObZGc+ye//uSJcnInJMCpwC6OQSJ/PCVS+wssdK2tZpywtfv0CtpFq8ueBNDQwqIvmWAqcAq1a2GqNuGsWse2dx1qln0WlCJ+KGxvHlz19GujQRkSMocKJA46qN+a7Td4y6cRTb9m6j6YimtH2vLanbUiNdmojIHxQ4UcLMuO2vt7Gi+wp6X9Wb6T9Np+7Aujw89WG2/7490uWJiChwok2pYqV4qslTpPRMoUODDvT/vj81X63JwDkDOZh+MNLliUghpsCJUpVPqsyw1sOY33k+9U+rT48pPag/uD5TUqZEujQRKaQUOFHu/Mrn83mHzxl36zgOpB+g1ahWtHynJUs3LY10aSJSyChwCgEzo02dNizttpR+Lfoxe+1s6g+pT9eJXdmctjnS5YlIIaHAKUSKxxSn1yW9SO2ZSre4bgydP5SaSTV55dtX2HdwX6TLE5Eop8AphMqXLk9SqySWdFtCk2pNeGz6Y9QdVJcPl32ogUFFJGQUOIVYnQp1mHj7RKbdOY3SxUpz89ibuXLElcxbNy/SpYlIFFLgCM3Pac6CzgsYct0Qlm9ezkVDL+LucXezdufaSJcmIlFEgSMAFC1SlM5xnUnpmcKjlz7K6CWjqTWgFolfJrLnwJ5IlyciUUCBI39StmRZ+jTvw/Luy2kV24qEmQnUSqrF24veJt3TI12eiBRg+S5wzCzBzNaa2cLg1CrTsifNLNXMVprZNZnmX2hmPwSXvWpmFpnqo8fZ5c5mbLuxfH3P15x+8ul0GNeBxsMa880v30S6NBEpoPJd4AT1d/fzg9NkADOrC7QHzgNaAoPMLCa4/mDgfiA2OLWMQM1R6fJqlzP7vtmMvGEk63ato8mbTbhl7C2s/m11pEsTkQImvwZOVtoAY9x9n7uvBlKBRmZ2OnCKu8/yQJ/ekcANEawz6hSxItzV4C5W9lhJwhUJTEqZRJ2BdXhi+hPs3Lcz0uWJSAGRXwOnh5ktNrM3zKxccF4V4NdM66wJzqsS/Hz4/COY2f1mlmxmyZs36wn7nCpTvAzxV8azqscq2tdrT59v+xCbFMvr817Xi99E5LgiEjhmNt3MlmQxtSFweewc4HxgPdA3Y7MsvsqPMf/Ime6vu3ucu8dVrFjxxHekkKpyShVG3DCCuf+YS63yteg8sTMNX2vIZz9+FunSRCQfi0jguHszd6+XxTTe3Te6+yF3TweGAo2Cm60Bzsz0NVWBdcH5VbOYLyEWd0YcX939FWPbjWX3/t20eKcF14+6nhVbVkS6NBHJh/LdJbXgPZkMbYElwc8TgPZmVsLMahDoHDDH3dcDu8yscbB3WgdgfFiLLsTMjJvr3syy7st4udnLfP3L1/x18F95YMoDbN2zNdLliUg+ku8CB3g52MV5MdAU6AXg7kuB94FlwKdAd3fPuHHQFRhGoCPBj4Be+hJmJYuW5NHLHiWlZwr3NbyPgXMHEpsUy3++/w/7D+2PdHkikg9YYR2sMS4uzpOTkyNdRtRasmkJD097mGk/TiP2L7G80vwVWtdujR6REinYzGyeu8flZtv8eIYjUaBepXp8esenTLp9EjFFYrjhvRto9nYzFm1YFOnSRCRCFDgSMmZGq9hWLO6ymKRrk1i4YSENX2vIfRPuY8PuDZEuT0TCTIEjIVcsphg9GvUgtWcqDzZ+kJGLRhKbFMuLX7/I3gN7I12eiISJAkfCplypcvS7ph9Luy3l6hpX8/SMp6kzsA5jlozRi99ECgEFjoRdbPlYxrUfx4wOMyhXshy3fXgbl71xGbPXzI50aSISQgociZimNZoy7/55DG89nNXbV9N4eGPu+OgOftnxS6RLE5EQUOBIRMUUiaFTw06s6rGKp5s8zUfLP6L2gNo8M+MZdu/fHenyRCQPKXAkXzi5xMm8cNULrOyxkrZ12tL7697EJsXyxoI3NDCoSJRQ4Ei+Uq1sNUbdNIpZ986i+qnVuXfCvcQNjePLn7+MdGkicoIUOJIvNa7amO86fceoG0exbe82mo5oStv32pK6LTXSpYlILilwJN8yM277622s6L6C3lf1ZvpP06k7sC4PT32Y7b9vj3R5IpJDChzJ90oVK8VTTZ4ipWcKHRp0oP/3/an5ak0GzhnIwfSDkS5PRLIpx4FjiVbGEi0mFMWIHEvlkyozrPUw5neeT/3T6tNjSg/qD67P5JTJenBUpAA47mjRlmhFgPbAHcBFwD6gBLAZmAy87vGeEuI685xGiy7Y3J0JKyfwyGePkLotlRbntKBvi77Uq1Qv0qWJRLVQjxb9BYFXPj8JVPZ4P9PjvRLQBPgeeMkS7c7cNC6SW2ZGmzptWNptKf1a9GPO2jk0GNKArhO7sjltc6TLE5EsZOcMp5jH+4ETXSe/0RlOdNm6ZysJXyYwOHkwZYqX4Zkmz/DAxQ9QomiJSJcmElVO5Awn2y9gs0Q7B1jj8b7PEu1KoD4w0uN9e24ajjQFTnRasWUFj0x7hEkpkzi73Nm83Oxlbjz3Rr34TSSPhOsFbB8ChyzRagLDgRrAqNw0KhIqdSrUYeLtE5l25zRKFyvNzWNv5oq3rmDeunmRLk2k0MtJ4KR7vB8E2gL/8XjvBZwemrJETkzzc5qzoPMChlw3hBVbVhA3NI6O4zqydufaSJcmUmjlJHAOWKLdBnQEJgbnFcv7kkTyRtEiRekc15mUnik8duljjFkyhloDapH4ZSJp+9MiXZ5IoZOTwLkHuATo7fG+2hKtBvBOaMoSyTtlS5alT/M+LO++nFaxrUiYmUDtAbV5e9HbpHt6pMsTKTSy3Wkg2qjTQOH1zS/f0GtqL5LXJRN3Rhz9r+nP5dUuj3RZIgVCuHqpXQYkAGcBRQED3OP97Nw0HGkKnMIt3dN5d/G7PPn5k6zdtZZ2ddvRp1kfapSrEenSRPK1cPVSGw70Ay4nMOJAXPCnSIFTxIpwV4O7WNljJQlXJDApZRJ1BtbhielPsHPfzkiXJxKVcnKGM9vj/eIQ1xM2OsORzNbuXMtTM55i5KKRVCxdkeebPs+9F9xL0SJFI12aSL4SrktqLwExwEcExlMDwON9fm4ajjQFjmQleV0yvab24ptfvqFepXr0a9GP5uc0j3RZIvlGuC6pXUzgMtqLQN/g9O/cNCqSX8WdEcdXd3/F2HZjSdufRot3WnD9qOtZsWVFpEsTKfDUS03kKH4/+DtJs5N44esXSNufRte4riRcmUD50uUjXZpIxITlkhqAJdp1wHlAyYx5Hu//zE3DkabAkezalLaJ+C/ieX3+65xS4hSe+9tzdG/UneIxxSNdmkjYheWSmiXaEOBWoCeBLtHtCHSRFolqlcpUYvD1g1nUZRGNqjTioWkPUW9QPcavGK8Xv4nkQE7u4Vzq8d4B+M3jPZHAqANnhqYskfynXqV6fHrHp0y6fRIxRWK44b0baPZ2MxZuWBjp0kQKhJwEzt7gzz2WaGcABwiMGC1SaJgZrWJbsbjLYpKuTWLhhoVc8NoF3DfhPjbs3hDp8kTytZwEzkRLtFOBV4D5wM/AmBDUJJLvFYspRo9GPUjtmcqDjR9k5KKRxCbF8uLXL7L3wN7jf4FIIZSrXmqWaCWAkh7vO/K+pPBQpwHJSylbU3j0s0cZv3I81cpWo0+zPtx63q168ZtEnZD2UrNEu8rjfYYl2o1ZLfd4/yg3DUeaAkdC4YvVX/DQtIdYuGEhjas2pv81/WlctXGkyxLJM6HupXZF8Offs5iuz02jItGqaY2mJP8jmeGth/Pz9p+5ZPgl3P7h7fyy45dIlyYScXrwUyREdu3bRZ9v+9B3Vl8AHr7kYZ64/AlOKn5ShCsTyb1QX1J76FjLPd775bhRs3YEXnVwLtDI3ZMzLXsSuBc4BDzg7lOD8y8E3gJKAZOB/3N3N7MSwEjgQmArcKu7/3y8GhQ4Ei6/7PiFJ6Y/weglo6l8UmV6X9Wbjg06ElMkJtKlieRYqC+pnRyc4oCuQJXg1AWom5tGgSXAjcBXmWeaWV2gPYHRDFoCg8ws47/KwcD9QGxwahmcfy/wm7vXBPoDfXJZk0hIVCtbjVE3jWLWvbOofmp17p1wL3FD4/jy5y8jXZpIWB03cDzeE4MPelYALvB4f9jj/WECZxRVc9Oouy9395VZLGoDjHH3fe6+GkgFGpnZ6cAp7j7LA6dkI4EbMm0zIvj5A+BqU9cgyYcaV23Md52+Y/RNo9m2dxtNRzSl7XttSd2WGtY6Jq2axDe/fBPWNkUgZ8/hVAP2Z/p9P1A9T6sJnDn9mun3Nfz/M6o1Wcz/0zbufhDYAWQ5uqKZ3W9myWaWvHnz5jwuXeT4zIz29dqzovsKel/Vm+k/TafuwLo8NPUhftv7W8jbd3fajW1HkzebcNP7N/Hjth9D3qZIhpwEztvAHEu0BEu0eGA2gTONLJnZdDNbksXU5hhtZHVm4seYf6xtjpzp/rq7x7l7XMWKFY9RhkholSpWiqeaPEVKzxQ6NOjAf77/D7FJsQyYM4ADhw6ErN1d+3ex9+BeGldtzNTUqdQdVJdHpz3K9t+3h6xNkQzZDhyP995AJ+A3YDtwj8f7i0dd372Zu9fLYhp/jGbW8Ofx2aoC64Lzq2Yx/0/bmFlRoCywLbv7JRJJlU+qzLDWw5jfeT71T6tPzyk9qT+kPpNTJodkYNCNuzcC0C2uG6t6ruKOv95B31l9iU2KZdDcQRxMP5jnbR7NxFUTeeGrF/RK70IkJ2c4eLzPA0YDHwNbLdGq5XE9E4D2ZlbCzGoQ6Bwwx93XA7vMrHHw/kwHYHymbToGP98MzPDC2tdbCqzzK5/P5x0+Z9yt4ziYfpDrRl1Hy3dbsmTTkjxtZ1PaJiAwAvYZJ5/BG23eIPn+ZM6reB7dJ3enwZAGfJr6aZ62eTR9Z/Xl2S+epearNXkt+bWwhp1ERk5eT9DaEi0FWA3MDP6ckptGzaytma0hMOL0JDObCuDuS4H3gWXAp0B3dz8U3KwrMIxAR4IfM7U9HChvZqnAQ8ATualJJNLMjDZ12rC021L6tejHnLVzaDCkAV0ndv0jKE5U5sDJcMHpF/BFxy/46JaP2HdwH9e+ey3XvnstyzYvy5M2j1VL3Blx1KlQhy6TutDwtYZ89uNnIW2zMFu2eRnrdq07/oohlJMznOeBxsAqj/caQDPg29w06u4fu3tVdy/h7qe5+zWZlvV293Pcvba7T8k0Pzl4Se4cd++RcRbj7r+7ezt3r+nujdz9p9zUJJJfFI8pTq9LepHaM5XuF3Vn6PyhxCbF8vK3L7Pv4L4T+u6sAgcCYdf23LYs7baUfzf/N7N+nUX9wfXpPqk7W/ZsOaE2j1VL3OlxzLx7Jh+0+4A9B/bQ4p0WXDfqukLxSu90T2fO2jkcSj90/JXzwN/e/Bs1X63Jc188x+79u8PS5uFyEjgHPN63AkUs0Yp4vH8BnB+askSkfOnyvHrtqyzptoQm1Zrw+PTHOXfguXyw7INc39/ZmBa4h1OxTNadZkoULcHDlz5M6gOpdInrwmvzXqPmqzXp+13fEw67zA6mH2Trnq1UKlMJM+OmujexrNsyXm72Mt/88g31BtWj5+SebN2zNc/aPJbd+3dzRt8zaD26Nau2rgpLm5+s/ISLh10clmey0vansXVv4N/7+a+ep1ZSLd5a+Bbpnh7Sdg+Xk8DZbol2EoGHNd+1RPsvoIuuIiFWp0IdJt4+kWl3TqNM8TK0G9uOK966guR1OR8pY1PaJsqVLHfc12NXKF2BAa0GsLjrYi4981Ie+ewRzht0Hh8v/zhPOjNs2bMFxzntpNP+mFeiaAkevexRUnumcv+F9zMoeRA1k2rSf1Z/9h/af4xvO3H/2/4/1u9ezyerPuG8Qefx4KcPsm1vaPse/bz9ZyBwTEL9TFbGmW3ClQl82+lbzix7JveMv4eLhl7EV//76jhb552cBE4bYA/Qi8D9lR8JDOApImHQ/JzmLOi8gCHXDWHFlhVcNPQiOo7ryNqda7P9HZvSNh1xOe1Y6lasy+Q7JvPpHZ9SomgJbnz/RpqOaMqC9Qtyswt/qgOOvLQHgbOvQdcNYnGXxVxc5eKwvNI7o54xN42h0/mdSJqTRM1Xa/Lq7FdD1k19U9omYiyGlJ4pIX8mK/O/96VnXsqse2fxTtt32JS2iSveuiJ8z2S5+3EnEoghgenZWbegTBdeeKGLFFTb9273x6Y95sWfL+6le5f2hC8SfPe+3cfd7oo3r/AmbzTJVZsHDh3wQXMGeYWXK7glmN8z7h5ft3Ndrr7rsx8/cxLwmT/PPO66k1dN9nMHnOsk4E3fauoL1i/IVZvHMvqH0U4CvnTTUnd3X7RhkTcb2cxJwGsl1fIJKyZ4enp6nrZ53/j7vPK/K//x+/pd6/2+8fe5JZiX71Pek2Yn+f6D+/OkrQkrJjgJ+Ny1c/80P21/mj8/83kv07uMF3++uD8y9RH/be9vx/wuINlz+Xc3W2c4Hu+HCLxaumxo409EsqNsybL0ad6HFd1XcF3sdSTMTKD2gNq8vejtY16X35i28U+XsXKiaJGidL2oKyk9U3j4kod5Z/E7xCbF8vzM59lzYE+OvutYZziHuzb2WhZ1WcSAaweweOPikLzS+/B66p9Wn2l3TuOT2z7BMFqPaU3zt5uzeOPivGtzzyZOK/P/j0XlkyoztPVQFnReQIPKDfL0mayMe3eZ2wMoXaw0z/ztmbA9k5WTS2q/Az9Yog23RHs1Y8rzikQk22qUq8H77d7n63u+5vSTT6fDuA5cPOzio46VtiltE5VKZ/+SWlZOLXkqr7R4heXdl3NNzWt47svnqD2gNu8ufjfbN6EzHkA9/A/g0RSLKUb3Rt1J6ZlCr8a98vyV3pvSNlHEivCXUn/5Y56ZcX2t6/mh6w/8t+V/mb9+Pg1fa8j9n9z/R/0n2mZWgdugcgOm3zWd8e3Hcyj9UJ48k5URqEfrLBKuZ7JyEjiTgGcJdBqYl2kSkQi7vNrlzL5vNiNvGMn6Xetp8mYT2o1tx+rfVv+xzoFDB9i2d1uO7uEcyzl/OYcPb/mQmXfPpFKZStz58Z1cMvwSvvv1u+NuuyltE0WLFOXUkqfmqM1ypcrR95q+LOu+jOZnN+fpGU9TZ2AdRv8w+oTOAjalbaJi6YoUsSP/JBaLKcYDFz9A6gOpPNDoAd5c+CaxSbG89M1L/H7w91y3uXH3xqMeCzOjde3WLOm2hP7X9P/jmawuE7vk6pmsTWmbOKXEKZQsWvKY64X8mazcXosr6JPu4Ui02r1vtyd8keCle5f24s8X98emPeY7ft/h63aucxLwQXMG5Xmbh9IP+VsL3vIz+p7hJOC3jr3Vf/7t56Ou32lcJz+j7xkn3O4Xq7/whkMaOgl442GNfdavs3L1Pa1Ht/b6g+tna92VW1Z669GtnQT8rP5n+ZgfxuTq/k6Z3mW816e9srXulrQt3nNyT49JjPFT/nWK9/mmj/9+4Pdst3XbB7d5zVdr5qi+fQf3ed/v+nrZf5X1mMQY7zaxm29O2xzaeziWaJ9Yov3dEq1YFsvOtkT7pyVap7yJPxE5UWWKlyH+ynhW9VhF+3rtefm7lwPP0gTfPJrbezjHUsSK0PH8jqzqsYrn/vYcE1ZOoPaA2jz1+VNZjpW2aU/OessdzZXVr2TuP+byRus3TuiV3jnpvVerfC3Gtx/P9LumU7ZkWdp/2J7L37ycOWvnZLu9tP1ppB1Iy/YlxRN9Jmtj2sZst5WheExxHrrkoSOeyToR2bmk9g+gCbDCEm2uJdpkS7QZlmirgdeAeR7vb5xQFSKS56qcUoURN4xg7j/mUrtC7T8CJ68uqWWlTPEyJDZNZGWPlbQ7rx3/+uZfxCbFMnTe0D89Ub9xd87/AB5NTJEY7ml4Dyk9U3imyTN8vOJjag+ozTMznsn2E/U57S4OcPXZVzP//vkM/ftQftz2IxcPu5g7P7qTX3f8etxtN+8JvB4lp23m9pms3OxfhsOfyTohOTkdIoHqJHAJCZxPAqVze1qVHyZdUpPCJD093T9Y+oHfPe5uT9ufFrZ2Z6+Z7ZcNv8xJwOsPru/Tf5zu7u5n9T/L7/rorpC0+b/t//PbP7zdScAr/7uyD58/3A8eOnjMbU568SR/cMqDuW5z5+87/cnpT3qJ50t4qRdK+bMznvVd+3Yddf3Za2Y7CfjElRNz3eaBQwd8yNwhXvHlik4C3uHjDr5mx5os1630SiXv/EnnXLeVGaHuFp05n4AyHu8LAbdEO/nE4k5EwiFj+Jg327xJ6WKlw9ZuoyqN+Pqer3nv5vfYuW8nzd5uRuvRrdmwe0PIzrSqla3Guze+y/f3fk+NU2sc95Xeew7sYff+3Sd0qfHkEifz4tUvsrLHStrUaXPc4WMyermdyL9B0SJF6RzXmdQHUnn8sscZs2QMtQbUIvHLRNL2p/2x3qH0Q2zZsyWkZ7bZlZPRov9B4BXOrwVnVQXGhaAmEYkiZsYt593C8u7Leenql/jy5y/Zd2hfyP8AXlz1Yr7t9C1jbhrDb3t/O+rwMTl5Juh4zjr1LEbfNJrvOn1HtbLVjjp8TEabeXE/7ZQSp/BSs5eOeCZr5KKRpHs6W/duJd3T8+wS5onIyRlOd+AyYCeAx3sKEPnIFJECoWTRkjx++eOk9Ezhpatf4q76d4W8TTPj1nq3srz7cl686sUsh4/Jy8DJcMmZl/Ddvd/x7o3vsjlt8xHDx/zxXEzpvHvz8OHPZHUc15GLh13Mx8s/BkJ77y67chI4+zze/xhBzxKtKEd5lbOIyNGcdtJpPH7545x+8ulha7NUsVI82eRJUnqm0LFBxz+90jtjLLq8/oNcxIpw+19vZ0WPFTzf9Pk/vdI7ZVsKJxc/mVLFSuVpm3DkM1ldJnUB8kfgmGejSx2AJdrLBF4t3QHoCXQDlnm8Px2y6kIoLi7Ok5NzPtquiBR8izYs4qFpDzFj9QzKFCtD2oE0fv6/nznr1LNC1ub6Xet5ZsYzvLnwTRznnHLnkPpAaEaHzpC2P42+s/oycdVEPr3z0z+NpJBbZjbP3eNytW0OAqcIcC/QAjBgKjDM47P5BfmMAkekcHN3Pln1CY9Me4SNaRvZ+MjG4z6JnxcWrF/Ak58/Se3ytfnvtf8NeXt5LSyB86eNEu0vQFWP97wbyS7MFDgiAoEhf3bs20GF0hUiXUqBcCKBUzTbjSTal0Dr4DYLgc2WaDM93h/KTcMiIvlBsZhiCpswyUmngbIe7zuBG4E3Pd4vBJqFpiwREYk2OQmcopZopwO3ABNDVI+IiESpnATOPwl0FEj1eJ9riXY2kBKaskREJNrkqtNANFCnARGRnAtXp4GSBLpFnwf80XfQ412vJhARkePKySW1t4HKwDXATAJjqe0KRVEiIhJ9chI4NT3enwXSPN5HANcBfw1NWSIiEm1yEjgHgj+3W6LVA8oC1fO8IhERiUrZvocDvG6JVg54FpgAnAQ8F5KqREQk6mQ7cDzehwU/zgTODk05IiISrXLSS60EcBOBy2h/bOfx/s+8L0tERKJNTi6pjQd2APOAfaEpR0REolVOAqeqx3vLkFUiIiJRLSe91L6zRFM3aBERyZXjnuFYov1A4FXSRYF7LNF+InBJzQD3eK8f2hJFRCQaZOeS2vUhr0JERKJedgJnI9AFqAn8AAz3eD8Y0qpERCTqZOcezgggjkDYXAv0DWlFIiISlbITOHU93u/0eH8NuBlocqKNmlk7M1tqZulmFpdpfnUz22tmC4PTkEzLLjSzH8ws1cxeNTMLzi9hZu8F5882s+onWp+IiOS97AROxhhq5OGltCUEXlX9VRbLfnT384NTl0zzBwP3A7HBKaOL9r3Ab+5eE+gP9MmjGkVEJA9l5x5OA0u0ncHPBpQK/p7RS+2UnDbq7ssBgicpx2VmpwOnuPus4O8jgRuAKUAbICG46gfAADMzL6xvlhMRyaeOGzge7zHhKCSTGma2ANgJPOPuXwNVgDWZ1lkTnEfw568A7n7QzHYA5YEth3+xmd1P4CyJatWqhWwHRETkSDkZaSBHzGw6gRe2He5pdx9/lM3WA9XcfauZXQiMM7PzCJxNHS7jDOZYy/480/114HUIvGL6WPWLiEjeClnguHuzXGyzj+A4be4+z8x+BGoROKOpmmnVqsC64Oc1wJnAGjMrSuA9PdtOoHQREQmBnAxtE3JmVtHMYoKfzybQOeAnd18P7DKzxsHeaR0IDCYKgXfzdAx+vhmYofs3IiL5T0QCx8zamtka4BJgkplNDS76G7DYzBYR6ADQxd0zzla6AsOAVOBHAh0GAIYD5c0sFXgIeCJMuyEiIjlghfVkIC4uzpOTkyNdhohIgWJm89w97vhrHilfXVITEZHopcAREZGwUOCIiEhYKHBERCQsFDgiIhIWChwREQkLBY6IiISFAkdERMJCgSMiImGhwBERkbBQ4IiISFgocEREJCwUOCIiEhYKHBERCQsFjoiIhIUCR0REwkKBIyIiYaHAERGRsFDgiIhIWChwREQkLBQ4IiISFgocEREJCwWOiIiEhQJHRETCQoEjIiJhocAREZGwUOCIiEhYKHBERCQsFDgiIhIWChwREQkLBY6IiISFAkdERMJCgSMiImGhwBERkbBQ4IiISFgocEREJCwUOCIiEhYRCRwze8XMVpjZYjP72MxOzbTsSTNLNbOVZnZNpvkXmtkPwWWvmpkF55cws/eC82ebWfXw75GIiBxPpM5wPgPquXt9YBXwJICZ1QXaA+cBLYFBZhYT3GYwcD8QG5xaBuffC/zm7jWB/kCfcO2EiIhkX0QCx92nufvB4K/fA1WDn9sAY9x9n7uvBlKBRmZ2OnCKu89ydwdGAjdk2mZE8PMHwNUZZz8iIpJ/5Id7OJ2AKcHPVYBfMy1bE5xXJfj58Pl/2iYYYjuA8lk1ZGb3m1mymSVv3rw5z3ZARESOr2iovtjMpgOVs1j0tLuPD67zNHAQeDdjsyzW92PMP9Y2R850fx14HSAuLi7LdUREJDRCFjju3uxYy82sI3A9cHXwMhkEzlzOzLRaVWBdcH7VLOZn3maNmRUFygLbTngHREQkT0Wql1pL4HGgtbvvybRoAtA+2POsBoHOAXPcfT2wy8waB+/PdADGZ9qmY/DzzcCMTAEmIiL5RMjOcI5jAFAC+Cx4f/97d+/i7kvN7H1gGYFLbd3d/VBwm67AW0ApAvd8Mu77DAfeNrNUAmc27cO2FyIikm1WWE8G4uLiPDk5OdJliIgUKGY2z93jcrNtfuilJiIihYACR0REwkKBIyIiYaHAERGRsFDgiIhIWChwREQkLBQ4IiISFgocEREJCwWOiIiEhQJHRETCQoEjIiJhocAREZGwUOCIiEhYKHBERCQsFDgiIhIWChwREQkLBY6IiISFAkdERMKi0L5i2sx2ASsjXUcIVQC2RLqIEIrm/YvmfQPtX0FX291Pzs2GRfO6kgJkZW7fy10QmFmy9q9giuZ9A+1fQWdmybndVpfUREQkLBQ4IiISFoU5cF6PdAEhpv0ruKJ530D7V9Dlev8KbacBEREJr8J8hiMiImGkwBERkbCI+sAxs5ZmttLMUs3siSyWm5m9Gly+2MwuiESduZGNfbvSzHaY2cLg9Fwk6swtM3vDzDaZ2ZKjLC+wxw6ytX8F9viZ2Zlm9oWZLTezpWb2f1msU2CPXzb3ryAfv5JmNsfMFgX3LzGLdXJ+/Nw9aicgBvgROBsoDiwC6h62TitgCmBAY2B2pOvOw327EpgY6VpPYB//BlwALDnK8gJ57HKwfwX2+AGnAxcEP58MrIqW//ZysH8F+fgZcFLwczFgNtD4RI9ftJ/hNAJS3f0nd98PjAHaHLZOG2CkB3wPnGpmp4e70FzIzr4VaO7+FbDtGKsU1GMHZGv/Cix3X+/u84OfdwHLgSqHrVZgj18296/ACh6T3cFfiwWnw3uY5fj4RXvgVAF+zfT7Go78H0V21smPslv3JcHT4ilmdl54SgubgnrscqLAHz8zqw40JPD/kjOLiuN3jP2DAnz8zCzGzBYCm4DP3P2Ej1+0D21jWcw7PKWzs05+lJ265wNnuftuM2sFjANiQ11YGBXUY5ddBf74mdlJwIfAg+6+8/DFWWxSoI7fcfavQB8/dz8EnG9mpwIfm1k9d898vzHHxy/az3DWAGdm+r0qsC4X6+RHx63b3XdmnBa7+2SgmJlVCF+JIVdQj122FPTjZ2bFCPwxftfdP8pilQJ9/I63fwX9+GVw9+3Al0DLwxbl+PhFe+DMBWLNrIaZFQfaAxMOW2cC0CHY46IxsMPd14e70Fw47r6ZWWUzs+DnRgSO99awVxo6BfXYZUtBPn7BuocDy92931FWK7DHLzv7V8CPX8XgmQ1mVgpoBqw4bLUcH7+ovqTm7gfNrAcwlUCvrjfcfamZdQkuHwJMJtDbIhXYA9wTqXpzIpv7djPQ1cwOAnuB9h7sXlIQmNloAj19KpjZGiCewM3LAn3sMmRj/wry8bsMuAv4IXgfAOApoBpExfHLzv4V5ON3OjDCzGIIBOX77j7xRP92amgbEREJi2i/pCYiIvmEAkdERMJCgSMiImGhwBERkbBQ4IiISFgocEREJCwUOCIiEhb/D2dzutVx3+rTAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax4 = plt.subplots()\n",
    "ax4.plot(f, np.unwrap(angles)*180/np.pi, 'g')\n",
    "ax4.set_ylabel('Phase(radians)', color='g')\n",
    "ax4.set_xlim([0, 3])"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
